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Isabella Karle Honored 


announced on November 17 at a 
jton that Dr. Isabella Lugoski Karle is the 
> Bower Award for Science, which consists 
1 $250,000. Dr. Karle is the first women to 
JS annual award of the Franklin Institute 
n 1834 to promote scientific inquiry and 
The award will be presented April 7, 1994 


Dr. Karle, wt a senicr scientist specializing in crystal- 
yraphy at the Naval Research Laboratory, is honored for her 
neering work in developing electron and X-ray diffraction 
hniques that make immediate determination of the three-di- 
sional structure of molecules possible. Her inventive work 
dged the gap between the theoretical formulas for crystal 
icture determination and the practical methodology re- 


quired to achieve useful applications. This has facilitated 
research in chemistry, biology, and medicine. Dr. Karie’s own 
studies have included crystallographic analyses of photorear- 
rangement products, ionophores, peptides, toxins, antitoxins 
cardioaciive substances and natural analgesics 


A native of Detroit, Michigan, Dr. Karle won a scholarship 
to the University of Michigan where she earned B.S., M.S. and 
Ph.D degrees with a specialty in physical chemistry. While at 
Michigan, she married Dr. Jerome Karle; and they both came 
to NRL in 1946. Dr. Jerome Karle is Chief Scientist of the 
Laboratory for the Structure of Matter at NRL. His work in 
crystal structure analysis was recognized by the 1985 Nobel 
Prize in Chemistry. Dr. Isabella Karle often collaborates with 
her husband on research projects 
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Remarks At 

The Bower Award Reception 
For Isabella Karle 

Navy Museum, Washington Navy Yard 
November 18, 1993 


Left to 


- 
Nora Slatkin 


right: Dr. Fred Saalfeld, Deputy Chief of Naval Research; RADM Marc Pelaez, Chief of Naval Research; the Honorable 
Assistant Secretary of Navy (Research, Development and Acquisition); Dr. Timothy Coffey, NRL Director of 


Research; Madeleine Tawney; Dr. Jerome Karle; Dr. Jean Karle; Dr. lsabella Karle: the Honorable Anita K. Jones, Director, Defense 
Research and Engineering; the Honorable John H. Dalton, Secretary of the Navy; Dr. Larry Tise, Executive Director, Franklin 


Memorial; Dr. James L. Powe 


President, Franklin Institute; and CAPT Paul G. Gaffney II, NRL, Commanding Officer. 





CAPT Gaffney: The nation’s Naval Research Laboratory 
(NRL) is honored to assist the historic and renowned Frank- 
lin Institute and Franklin Memorial, this evening, in its 
formal announcement of the prestigious Bower Award and 
Prize for Scientific Achievement. 

We have a very short program featuring several top 
leaders in national science and defense policy. So, without 
further ado, it is my distinct pleasure to introduce to you, 
President and Chief Executive Officer of the Franklin Insti- 
tute, Dr. James L. Powell, and Dr. Larry Tise, the Executive 
Director of the Franklin Memorial. 

Dr. Powell: Thank you very much Captain Gaffney. Good 
evening everyone. 

It was a former Assistant Secretary of the Navy, Franklin 
Delano Roosevelt, who said “When you're giving a talk after 
dinner be sincere, be brief, be seated.” So, I’m going to try to 
follow his advice this evening, but it’s a great honor for me 
personally to be here as a scientist. I’ve known about the 
important work of NRL and its various incarnations all my 
career and I’m very privileged to be here. 
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The Franklin Institute was founded in 1824 to promote 
the “mechanic arts,” as they referred to them in those days, 
and we've been doing science, in one way or another, ever 
since. One of the things that those of us who are concerned 
about science realize is that we do not have a very favorable 
image. The media does not help. Instead of encouraging 
children to go into science, they seem to have just the 
opposite effect, and so I think we need to seize every oppor- 
tunity that we have to recognize scientific achievement. 

I was very pleased when I got to the Franklin Institute a 
few years ago to find out that we have a distinguished 
program for recognizing scientific achievement. And, when 
that achievement goes to a woman or a member ofa minority 
group, so much the better, because those are the people we 
especially need to encourage into science. In just a moment, 
Dr. Tise, my colleague, will say a bit about the process by 
which the awardee is selected, and you will see that it is a 
result of an international search whose ground rules are only 
that we try to find the very best person that we can in science. 
Just to show you the conclusion that we came to about Dr. 





Karle, before I turn to Dr. Tise, I want to read you part of the 
citation that will be presented to her and is the basis for this 
award. “The Bower Award and Prize in Science goes to 
Isabella L.Karle, physical chemist and crystallographer, for 
pioneering contributions in determining the three dimen- 
sional structure of molecules, making use of both X-ray and 
electron diffraction, and, in particular, for her definitive 
introduction of a procedure to reveal molecular structure 
directly from X-ray studies. This inventive work has pro- 
foundly facilitated studies in chemistry, biology, and medi- 
cine, permitting routine and immediate determinations of 
molecular structures. She has thrived in her profession and 
served as an inspiration to her colleagues and students with 
a curiosity and creativity that make her a kindred spirit of 
Benjamin Franklin himself.” So, we are very pleased to be 
able to make this award to Dr. Karle. And now, I'd like to turn 
to Larry Tise. 


The Honorable John H. Dalton, Secretary of the Navy. 


Dr. Tise: Not as a scientist but an historian, I was very 
pleased, when I came to the Franklin Institute some five 
years ago, and had the responsibility of establishing the 
procedures for the Bower Award, to find that my chief mis- 
sion was to find that individual in our world today, every year, 
who reflects the scientific and humanitarian genius of Ben- 
jamin Franklin. 

In setting up the process of finding this person every 
year, we established a worldwide nomination process, an 
international board of advisors; we established a selection 
process that is more rigorous, or equally rigorous, to any 
other award in the world and, every year, we receive thou- 
sands of inquiries about this award. After discouraging many 
people from submitting a nomination, we receive between 
50 and 100 fully documented nominations, which are then 
reviewed on three levels; locally, among scientists in our 
community and then internationally, among scientists who 
are expert in the field of the nomination, and then, finally, 
we bring together an international selection committee of 
people, mainly not from the United States, to help us decide 
who should receive the award in this particular year. 

This year Isabella Karle was selected not only because of 
an outstanding scientific achievement that you just heard Dr. 
Powell describe, but also because her work has applications in 
so many fields. I wish you could have heard the selection 
committee when they discussed this nomination and made the 
selection. They said, “Karle and Karle.” It is known to everybody 
all over the world and because of her life and contributions, not 
only as a scientist, but also as a citizen and, particularly in this 
year, as an inspiration to women to choose careers in science. 
It has been my pleasure to sort of guide and direct this process 
this year and to come to know Dr. Karle. 

I would like to invite all of you to Philadelphia on April 
7th and 8th, next year, when we will not only present her with 
a check and with a gold medal, but we will be holding two 
days, actually three days, of events in which we will be 
exploring the work that she has done and that others in her 
field have done, and we will be looking at related issues, such 
as the role of women in science today. We will celebrate her 
work and, of course, we'll spend some time talking about 
Benjamin Franklin as well. Thank you. 

CAPT Gaffney: Ladies and gentlemen, in his first ap- 
pearance in connection with science and technology, the 
Honorable John H. Dalton, Secretary of the Navy. 

Secretary Dalton: Thank you very much, and I want to 
welcome each of you here tonight to this beautiful museum. 
Several people, during the cocktail hour, mentioned to me 
that this is the first time they had been in this building and 
I hope that you'll come back and take it in. It really is a 
wonderful place and we’re fortunate to have it. 

I'm very pleased to be here tonight on this occasion. It’s 
a tremendous day for the Department of the Navy and it’s 
great that we’re here to honor Dr. Isabella Karle, as the 
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recipient of the 1993 Bower Award and Prize for Achievement 
in Science. The Bower Award, which recognizes outstanding 
achievement in the life or physical sciences, serves as a 
memorial to that great American scientist and humanitarian, 
diplomat and patriot, Benjamin Franklin. Dr. Karle’s scien- 
tific achievements stand solemnly within the tradition estab- 
lished by Benjamin Franklin. A tradition that combines 
scientific excellence with public service. The great Franklin 
once advised a friend to continue with his “excellent experi- 
ments...produce facts, improve science and do good to man- 
kind.” Dr. Isabella Karle has excelled in all these regards. An 
outstanding researcher, laboratory teacher, lifelong civil servant 
and parent, she has set an example to which her colleagues and 
succeeding generations of scientists can aspire. She has indeed 
carried out excellent experiments, produced facts, improved 
science, and done good to mankind. At the same time, she has 
raised a family of achievers in their own right. 

As Secretary of the Navy, I’m especially pleased that the 
Department of the Navy can claim Isabella as its own. Her 
work stands as an example for the high quality research 
carried out by the men and women at the Naval Research 
Laboratory. We, in the Navy, are fortunate to have developed 
and maintained a national resource like NRL, which can 
attract and retain world class scientists, like Isabeila and her 
husband, Nobel Laureate, Dr. Jerome Karle. Dr. Isabella 
Karle’s career spans over 47 years with the Navy, at the Naval 
Research Laboratory. After receiving her Ph.D. in chemistry, 
at the age of 22, she and her husband both worked on the 
Manhattan project. Both were subsequently appointed to 
NRL. Isabella’s efforts for our Navy and our nation have 
combined hard work with insights into the processes of 
nature. Dr. Karle’s master work is her now famous symbolic 
addition procedure for determining molecular structures, an 
advance that has greatly aided medical and pharmaceutical 
research. The author of almost 300 scientific papers and the 
recipient of numerous awards and citations, including the 
Secretary of the Navy’s Distinguished Achievement Award in 
Science, Dr. Isabella Karle is the only woman ever elected to 
the Chemistry Section of the National Academy of Sciences 
and now she is the first woman and first US government 
scientist to win the Bower Award and Prize. As the Depart- 
ment of the Navy right-sizes and reshapes its forces with a 
post-cold-war world, science and technology become ever- 
more important in continuously improving our capabilities 
to use the sea and preserve peace. 

Great institutions such as the Naval Research Labora- 
tory and the Franklin Institute as well as dedicated scientists 
like Dr. Karle are critical if we are to be successful in these 
efforts. They will provide the foundation for America’s future. 

On behalf of the President of the United States and the 
Department of the Navy, I am very proud to be here to 
personally congratulate Dr. Isabella Karle on her selection to 
receive the Franklin Institute’s Bower Award and Prize for 
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The Honorable Anita K. Jones, Director, Defense Research 
and Engineering 





Achievement in Science. Isabella, we’re proud of you. Con- 
gratulations and God bless you. 

CAPT Gaffney: Thank you Mr. Secretary. Ladies and 
gentlemen, the champion of defense science and technology, 
and I’m very pleased to say a rather frequent visitor to the 
Naval Research Laboratory. It’s my pleasure to introduce the 
Director of Defense Research and Engineering (DDR&E), the 
Honorable Anita K. Jones. 

Dr. Jones: Thank you Capt. Gaffney, and Secretary Dal- 
ton, Dr. Powell and Dr. Karle. 

It’s an honor to join with this group of guests to cele- 
brate the announcement of the award of the Bower Prize to 
Dr. Isabella Karle. It’s a particular pleasure to recognize this 
prize going to a renowned scientist from within the Depart- 
ment of Defense community. 

My first speech when I became the DDR&E was at the 
Naval Research Laboratory. It was at a seminar on contem- 
porary science honoring Dr. Jerry Karle on the occasion of 
his 75th birthday, and now I have the pleasure of being here. 
When I walked in, I had intended to ask whether there are 
some other members of the family that I should be on the 
lookout for. I was sitting next to Dr. Karle and she pointed 





Dr. Isabella Karle, Naval Research Laboratory Sen 
and recipient of the 1993 Bower Award 





out that the radiation-related material depicted in the left- 
hand panel of the display over here was the result of collabo- 
ration with a daughter. So | did not need to ask. 

An award like this, among other things, recognizes 
quality, and quality is crucial in a scientific enterprise. It’s 
particularly my job, anda fine one I may add, to try to defend, 
to catalyze, and to find new opportunities for injecting qual- 
ity into the science and technology program of the Depart- 
ment of Defense. Also it’s a real pleasure to see work of the 
highest distinction recognized by an award as eminent as the 
Bower Award from the Franklin Institute. Such an award 
spills over to the colleagues and to the environment in which 
a scientist works. I think the Naval Research Laboratory can 
be very proud of its support and the conditions it made 
possible, and which played a part in this award as well. 

It’s very difficult to express to someone who's not a 
scientist how exciting it can be and what a difference a single 
individual can make. Also it is difficult to make a description 
of Dr. Karle in down-home terms - I’m from Texas. One of 


the things she is noted for is creating electron and X-ray 
diffraction techniques, that allow one to see 3-D structures 
that you couldn’t see before and, so, a way of describing her 
work is: “she gave eyes to scientists across the world to see 
what they could not see before her.” To see into what was 
darkness before — that is a very exciting experience. Dr. Karle 
is the first woman, as Secretary Dalton said, to receive this 
Franklin Institute Award and it’s a particular pleasure for me 
to join this company, in recognizing her, and to know that 
we have such high recognition of a woman who is a role 
model for both men and women, but particularly women, 
who are moving through the defense laboratory science 
community. 

It is a pleasure to offer my congratulations to you Dr. 
Karle. I'd like to close by bringing to you a message from Les 
Aspin, the Secretary of Defense. He writes, 

“It is my pleasure to congratulate Dr. Isabella Karle for 
winning the 1993 Franklin Institute Bower Award and Prize 
for Achievement in Science. This award is a testament to Dr. 
Karle’s outstanding career as a physical chemist and crystal- 
lographer. Her work at the Naval Research Laboratory has 
enormous implications for the advancement of chemistry, 
biology, and medicine. Not only is Dr. Karle the first woman 
to win the highly prestigious Bower Science Award, she is 
also the first recipient from the community of Defense De- 
partment laboratories. We’re very proud of her. Defense 
research has been very successful in making our armed 
forces the world’s best. Dr. Karle’s recognition also under- 
scores the potential for private sector uses of the technolo- 
gies developed in our military laboratories. 

“In honoring Dr. Karle and her work, the Franklin 
Institute is highlighting the benefits of technology transfer 
from the military to the civilian sector. We are deeply grateful 
for that recognition by the institute. This ceremony provides 
a wonderful means of expressing gratitude to Dr. Karle. I 
appreciate the opportunity to add my thanks for her many 
remarkable contributions to science and to our nation.” 

CAPT Gaffney: Ladies and gentleman, reflecting on her 
magnificent career, addressing Navy and national scientific 
challenges, the senior scientist in NRL’s Laboratory for the 
Structure of Matter is the 1993 Bower Award selectee, Dr. 
Isabella Karle. 

Dr. Karle: Thank you Dr. Powell, Secretary Dalton, Dr. 
Tise, Secretary Jones for this magnificent event. I am indeed 
honored, excited, and grateful to have been named the re- 
cipient of the Bower Award. I wish to thank the trustees of 
the Franklin Institute and the selection committee for all 
their hard work. I also deeply appreciate having been nomi- 
nated by Dr. Coffey, the Director of Research, the Naval 
Research Laboratory and Dr. Rath, the Associate Director for 
Material Science. 

The Naval Research Laboratory deserves much credit for 
supporting basic research and for providing a suitable and 
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The Honorable Nora Slatkin, Assistant Secretary of the Navy 
(Research, Development and Acquisition) to the left and Lr 
Isabella Karle 





exciting environment for science to flourish. When my hus- 
band and I first arrived at the Laboratory 47 years ago, we 
were immediately thrust into an atmosphere of exciting 
research and surrounded by exceptional people, whose 
names still appear on the front page of the Laboratory news- 
paper and one of whom is here at this gathering. Dr. Herbert 
Friedman was our immediate supervisor and, as most of you 
know, he has won world acclaim for spatial astronomy and 
was presented with very many awards, amongst them the 
Wolf Prize of Israel of 1987. Dr. Richard Tousey was just a few 
doors down the hall and he, in 1990, received the Hale Prize 
for solar spectroscopy. Somewhat later, Dr. Timothy Coffey 
and Dr. Bhakta Rath, of course, they’re younger so they had 
to be later, came to the Laboratory and they too have pres- 
tigious awards. Dr. Coffey received the 1991 Franklin Insti- 
tute Delmer Fahrney Medal for leadership in science and 
technology. Dr. Bhakta Rath received the 1992 Minerals, 
Medals, and Materials Society Award for Materials Science. 
Almost anytime we open Labstracts, which is the Laboratory 
newspaper, there is someone who is receiving an important 
award. Last week, I think it was, Dr. Robert Lehmberg and 
Dr. Richard Obenschain received the American Physical So- 
ciety Award for Excellence in Plasma Physics. So, you see, 
these are only a very few examples of the national and 
international recognition accorded to NRL scientists, and 
that makes NRL an important and interesting place to work, 
especially since all of these people received their awards in 
rather diverse fields of science. The mood, prevalent in the 
Laboratory when we first came in the late 1940s, can be 
measured by the following occurrence. Jerome designed an 
electron diffraction apparatus, which was very large and 
complicated and I assembled the pieces, aligned and fine- 
tuned all the components, and eventually collected and ana- 
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lyzed data. Well, during the process of building the appara- 
tus, when the on-site foundry and machine shops were 
supposed to cast and fabricate the large brass segments, 
Captain Brad Bennett, (he must not have been a Captain then, 
he was a Commander, and at that time, he was in charge of 
Support Services) suggested and authorized that two of every- 
thing be fabricated. This was an excellent suggestion because it 
greatly facilitated the assembly of the apparatus. We had a 
successful first instrument much more quickly because if one 
piece didn’t work the other one did. Eventually, we had two 
apparatuses. So, as a result, we had papers with experimental 
data already analyzed, and in print, barely two years from the 
time we first started designing the instrument. 

A popular slogan in the 40s was “Join the Navy and See 
the World.” I didn’t know it applied to civilian scientists, but 
apparently it did because it turned out that way. 

Our subsequent work concerned X-ray diffraction and 
the elucidation of the three-dimensional structure of mole- 
cules. A key problem existed; in fact, half of the necessary 
information seemed to be lost in the process of the diffraction 
procedure. It was called the phase problem. If we could 
evaluate the phases as well as the intensities, then we could 
work backwards and determine the structures of these mole- 
cules. Well, Jerome Karle, with a colleague, Herbert 
Hauptman, who worked in our group at the time, solved the 
problem of obtaining the phases theoretically. The results 
were presented in terms of an infinite set of inequalities, 
probabilistic functions, and transcendental functions. It was 
I who bridged the gap between the theory and experiment by 
translating these mathematical results into practical proce- 
dure. Jerome always said he liked working with a paper and 
pencil much better than anything else. There was a great 
amount of trepidation when we gambled our research budget 
to have IBM do a final refinement on my first results. You 
see, at that time, this was the mid-1950s — IBM was the only 
company that had a computer that was capable of doing the 
problem that we needed to be done. Their computer was 
displayed in New York City on Fifth Avenue behind a plate- 
glass window, so that passerbys could see this wondrous 
machine doing its thing. Well, fortunately, their machine 
worked. My results worked and that was the beginning of 
doing several hundred (I’ve not counted them) structures 
over the years, in various areas of science. The successes 
we've had, in this sort of activity, have spawned numerous 
advance-study workshops that took place around the world 
for the next 20 to 25 years. These are the places I remember; 
there were a number of others. Cambridge University — | 
think that was the first one — and then Oxford, York, Parma, 
Erice (on top of a mountain in Sicily) Prague, Warsaw, 
Hungary, Kyoto, Brazil and so forth. In that way, we saw the 
world rather thoroughly. However, I didn’t get to the North 
Pole until two years ago. There were a number of scientists 
there, and I presented some lectures. I should say that this 





was in Tromso, at the university, way up above the Arctic 
Circle. The notes and instructions that I provided for the 
workshops became teaching materials in various universities 
and the basis for a number of computer programs. These 
courses were instrumental in expanding the number of crys- 
tal structures solved around the world, from a handful each 
year to much more than 10,000 a year, at present. We know 
how many because these structures are all in a data base 
called the Cambridge Data Base and they are having a diffi- 
cult time keeping up with the exponential number of pieces 
of data that they have to include every year. 

Personally, I’ve had a satisfying and fascinating career. 
Much of the success is due to sufficient support. At present, 
with our economy at a lower level, there is a tendency to 
de-emphasize basic research in many institutions. I would 


like to mention two recent publications that address the 
profits of society from basic research. The first is the Physics 
Today issue. In fact, it’s the October 1993 issue that is devoted 
to the benefits of basic research and some of the titles are: 
“From Basic Research to High Technology,” such research is 
the seed corn of the technological harvest that sustains 
modern society; “The Birth of the Laser Era.” and “Fiber 
Optics,” these are the components that went into making 
fiber optics practical. Fundamental research in glass science, 
optics, and quantum mechanics have matured into a tech- 
nology that is now driving our communications revolution. 
Then, of course, there are the magnetic resonance experiments, 
superconductivity, semiconductors, cyclotrons, and so on and 
so on. A somewhat similar publication came out recently from 
the American Chemical Society, under the direction of Ernest 


Dr. Isabella and Dr. Jerome Karle with their daughters Madeleine Tawney on the left and Dr. Jean Karle on the right at the ceremony 
Another daughter Louise Hanson was not able to attend 
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Eliel. It’s called Science and Serendipity, and it has a number 
of fascinating articles on how basic research, that was not 
directed toward anything, has resulted in nylon, teflon, and the 
platinum-plated bullet (renamed cisplatin, an anticancer drug). 
For the medical field, chemistry has provided cortisone, designer 
genes, AZT, and so on. I have mentioned chemistry and physics. 
Of course, similar examples exist in medicine and physiology, 
biology, agriculture and all the other scientific areas. Everywhere 
we look in modern society we see the end-product of basic 
research. Often it takes 10 or 50 or even more years for the 
knowledge derived from basic research to be applied to high 
technology, but the pool containing basic information must be 
constantly replenished, otherwise it will run dry. 

Again, | want to express my thanks to the Franklin 
Institute for choosing me as the Bower Award recipient and 
to the Navy and the Naval Research Laboratory for all of their 
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support for the last 47 years. Thank you. I also want to thank 
everybody for helping me celebrate tonight. 

CAPT Gaffney: Thank you Dr. Karle, Dr. Jones, Mr. 
Secretary, Dr. Tise, Dr. Powell and especially to the Franklin 
Institute, today. This is a memorable evening for the Navy, 
the Naval Research Laboratory and I’m sure for all the Karles 
that are here. We all look forward to the grand ceremony at 
the Franklin Memorial in Philadelphia in the spring. 

Now, if all of you will help us with the rest of the 
reception, we'll let the press spend some time with Isabella 
alone. 

Thank you very much everyone. 








Molecular Structure: Four Articles 


But beyond the bright search lights of science 
Out of sight of the windows of sense, 

Old riddles still bid us defiance, 

Old questions of Why and of Whence. 


When Science from Creation’s face 
Enchantment'’s veil withdraws, 
What lovely visions yield their place 
To cold material laws! 


Thomas Campbell (1777-1844) — To The Rainbow 


The following four articles by Isabella Karle represent three important and distinct 
areas of science which she developed and which have changed our thinking about the 
structure of matter. By developing X-ray diffraction techniques for determining the 
three-dimensional structure of molecules, she has made possible new avenues of research 
in chemistry, biology and medicine. 

The three areas of science and the accompanying articles are: 


Structure analysis of molecules in the vapor state by electron diffraction — 
“Internal Motion and Molecular Structure Studies by Electron Diffraction,” 
(1949). 


The determination of phases directly from the measured intensities of X-ray 
reflections — “Crystal and Molecular Structure of p,p’-Dimethoxybenzophenone 
by the Direct Probability Method” (1957) and “An Application of the Symbolic 
Addition Method to the Structure of L-Arginine Dihydrate,” (1964). 


Elucidation of molecular formulae and determinations of conformations of 
peptides, steroids, and alkaloids — “Crystal Structure of [Leu]zervamicin, a 
Membrane Ion-Channel Peptide: Implications for Gating Mechanisms,” (1991). 


In order to maintain authenticity and accuracy, the four articles have been photo- 
graphically duplicated as originally printed. 
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Reprinted from Tue JourNAL oF CuemicaL Puysics, Vol. 17, No. 11, 1052-1058, November, 1949 


Printed in U 


S. A. 


Internal Motion and Molecular Structure Studies by Electron Diffraction* 


ISABELLA LuGoski KARLE AND JEROME KARLE 
Naval Research Laboratory, Washington, D. C. 


(Received January 31, 1949) 


A procedure has been developed for the determination of molecular structure by electron diffraction 
which yields accurate intensity data and obviates the necessity for visual examination of the diffraction 
photographs. The theory for computing radial distribution curves has been extended to permit accurate 
curves to be obtained from scattering data covering only a restricted range of angle. From this method, it 
is possible to obtain not only equilibrium distances but also the probability distributions for the vibrational 
motion between pairs of atoms in a molecule. The procedure has been applied to CCl, and CO; and when 
comparisons may be made with spectroscopic results, satisfactory agreement is obtained. 


HE study of molecular structures using electron 

diffraction by gases has usually involved the 
visual examination of the diffraction pattern' which 
consists of small oscillations about a steeply falling 
background. The positions of the apparent maxima and 
minima are measured, their relative intensities are 
estimated visually, and some estimate is made of the 
shape of those features which appear to be neither 
maxima nor minima. The observed curve is used as the 
basis for evaluating the equilibrium distances in the 
molecule being studied. This is accomplished by calcu- 
lating a radial distribution curve which is used as a 
guide for selecting various models to represent the 
structure of the molecule. From these models intensity 
curves are computed and compared to the visually 
observed curve. The model chosen to represent the 
structure of the molecule depends upon the agreement 
between the computed and observed curves. Some 
differences are found in the form of the functions 
used by various investigators to compute the radial dis- 
tribution and intensity curves. 

The uncertainties in the visual examination of the 
diffraction patterns and in the mathematical analysis, 
make it very desirable to establish a quantitative pro- 
cedure which uses an unambiguous theory and elimi- 
nates the visual study of the photographs. A procedure 
to eliminate this visual study has been developed by 
the Norwegian school and they have used it to study 

* Presented at the meeting of the American Society for X-Ray 


and Electron Diffraction, Columbus, Ohio, December 16, 1948. 
'L. O. Brockway, Rev. Mod. Phys. 8, 231 (1936). 


10 Naval Research Reviews 


the structures of a large number of molecules.* Their 
procedure involves the use of a rotating sector’ in order 
to accentuate the oscillating features of the pattern and 
the photographs are scanned with a microphotometer. 
The main difference between their procedure and the 
one to be outlined concerns the method of analyzing 
the scattering data. By our method accurate values 
of the equilibrium distances can be expected and, in 
addition, it is possible to evaluate the magnitude of the 
vibrational motion between pairs of atoms in a molecule. 


THEORY 


When free molecules are struck by a beam of fast 
electrons, part of the scattering of the primary beam 
arises from the spacings between the atoms in indi- 
vidual molecules and can be represented by the formula 


Ials)@EE cud. (ivej) (1) 


i=l j=l 


The rest of the scattering does not depend upon the 
molecular structure and forms a steeply falling back- 
ground upon which the molecular scattering, Eq. (1), 
is superimposed. The upper limit » in the summation is 


2 (a) Chr. Finbak, Avhandl. Norske Vid.-Akad. Oslo, I Mat.- 
Naturv. Kl., No. 7 (1941). (b) Chr. Finbak and O. Hassel, Archiv. 
f. Mat. o. Naturvid. B45, No. 3 (1941). (c) H. Viervoll, Acta Chem. 
Scand. 1, 120 (1947). (d) O. Hassel and H. Viervoli, Acta Chem. 
Scand. 1, 149 (1947). 

*(a) C. Finbak, Avhandl. Norske Vid.-Akad. Oslo, Mat.- 
Naturv. KI., No. 13 (1937). (b) P. P. Debye, Physik. Zeits. 40, 66, 
404 (1939). 
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equal to the total number of atoms in the molecule. The 
coefficients cj; are characteristic of the ith and jth 
atoms and in the subsequent analysis are assumed to be 
constant. In practice, this is a good approximation, 
except for very small scattering angles, if the molecular 
scattering is obtained by dividing the total scattering 
by the background scattering. If a molecule is absolutely 
rigid, then the interference function A,; is equal to 
sinsr;;/sr;; where r;; is the distance between the nuclei 
of the ith and jth atoms and s= (4 sin@/2)/X where @ 
is the scattering angle and X is the wave-length of the 
electron beam. Since there is some motion of the nuclei 


in all molecules, the interference function A,; assumes 
the form 


A w= f P;;(p)(sinsp/sp)dp, (2) 
0 


where P;;(p)dp is the probability that the distance be- 
tween the ith and jth atoms has a value between p and 
p+dp. 

A Fourier transformation may be made on the 
molecular scattering function, 


=F E cy f Pij(p)(sinsp/sp)dp, (ij) (3) 
0 


i=1 j=1 
by multiplying it by ssinsr and integrating with re- 
spect to s. The result of such a procedure is a function 


of the variable r which represents the radial distribution 
curve, D(r) ; 


a 


b= f sI m(s) sinsrds (4) 
0 


nn “ ” Pii(p) 
=> > caf sinsrds f : 
i=1 j=l 0 0 wv 


7s 8 Pi;(r) 
res > Cir , (4b) 


2 i=l j=1 r 


sinspdp (4a) 


The ordinates of the radial distribution curve essen- 
tially represent the probability of the occurrence of the 
internuclear distances, r, in the molecule. It has been 
found in our experiments that in the vicinity of an equi- 
librium distance, 7;;, which is not very close in value 
to other ones in the molecule, the function, D(r) can 
usually be represented by 


D(r) = ¢;;(rh;;)} exp[ —hi(r—7:5)? )/ 20:5, (5) 


where (2h;;)-? is a measure of the average deviation 
from the equilibrium position and is equal to (I;7)+ 
where /,;; is the displacement from equilibrium. The 
quantity r in the denominator in Eq. (4b) may be 
replaced by 7;; when 4,; is sufficiently large. 


‘The derivation follow the treatment by P. Debye, J. Chem. 
Phys. 9, 55 (1941). 


DIFFRACTION STUDIES 


In order to evaluate the integral in Eq. (4), the molecu- 
lar scattering data must be known to values of s equal 
to infinity. In practice, only a limited amount of data are 
available. Nevertheless, this limited data uniquely deter- 
mines the structure of a molecule. The procedure used 
is to compute the following modified function,® 


so)= f » sI,(s) exp(—as*) sinsrds, (6) 
0 


where Smax bounds the range of scattering data and a 
can be so chosen that the integrand of Eq. (6) makes 
no essential contribution beyond Smax. The computed 
f(r) curve, therefore, is the same as that which would 
have been obtained if sna. approached infinity. It is 
possible then to obtain the desired function, D(r), from 
f(r) by making use of the relationship* 


1 ) 
fms f D(o) expl—(r—p)*/4a dp, (7) 


2(xa)i/_.. 


where a is the same as that used in Eq. (6). The intro- 
duction of D(r) as defined in Eq. (5) leads to an evalua- 
tion of the integral in Eq. (7), giving for f(r) the ex- 
pression 


f(r) - ci;(wh;;)3 
Xexpl — his(r—1:5)*/ (4ahis+ 1) ]/27:;(40h;;+1)'. (8) 


In practice, the bell-shaped peaks which result from 
the computation of Eq. (6) are curve fitted by the func- 
tion, B;; exp — Hi,(r—1;;)"], where Hi;=hi;/(4ah;;+ 1) 
(see the exponent in Eq. (8)). Since a is the value used in 
computing Eq. (6), 4;; can be readily obtained. The 
values of h;; so determined, as pointed out, are the 
measure of the magnitude of the vibrational motion 
between pairs of atoms in a molecule. 

Equation (§) may be used to illustrate the unique 
determination of the function, D(r), from a limited 
amount of data. By this we mean in a mathematical 
sense that the data may be extended uniquely beyond 
its experimentally determined range. There is, of course, 
always some uncertainty in experimental data. In a 
practical sense, then, it is possible to extend a function 
on the basis of a uniqueness property in a definite 
fashion beyond its known experimental range, but the 
uncertainties in the experiment will introduce some 
uncertainty into the extension. There is a unique rela- 
tion between the function in Eq. (8) whose parameters 
may be determined from experiment and the desired 
function in Eq. (5). If 4;; and 7;; are evaluated by ap- 
plying Eq. (8) to a computed curve, then D(r) is a 
known numerical function of r. From one rather re- 


* The damping factor exp(—as*) was introduced by Degard in 
order to improve radial distribution calculations. See C. Degard, 
Bull. Soc. Roy. Sci. Liege 12, 383 (1937). The effect of this and 
other damping factors on radial distribution curves has been 
studied by J. Waser, Thesis, California Institute of Technology. 

SE. C. Titchmarsh, Introduction to the Theory of the Fourier 
Integral (Clarendon Press, Oxford, 1937), p. 51. 
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Fic. 1. Typical microphotometer traces of a light and medium 
exposure of CCl,. The vertical scales are not the same. 


stricted, but simple, point of view, this uniqueness can 
be based upon the theoretical result’ that the computed 
radial distribution curve can be represented by the 
sum of Gaussian distributions. From another point of 
view, the study of the mathematics of positive Fourier 
integrals shows that the scattering function, /(s), known 
only in a restricted range of angle is related within 
narrow bounds to only one positive distribution func- 
tion. This idea has been used in our laboratory to im- 
prove the accuracy of determining electron distribu- 
tions*® about atoms. The restrictiveness of positive 
functions has also been found useful for small angle 
scattering and in the evaluation of the phases of the 
Fourier coefficients in the determination of crystal 
structure.* 

In the case of gas molecules, though, it is not neces- 
sary to resort to powerful mathematical arguments to 
obtain accurate information from a restricted amount of 
data. If the experimental curve, f(r), is merely fit well by 
a function in the vicinity of 7;; (though the function may 
even become negative at some larger distance), a good 
evaluation of the desired D(r) in the vicinity of 7;; can 
be obtained. This may be seen from Eq. (7) which 


?R. W. James, Physik. Zeits. 33, 737 (1932). 


*(a) H. Hauptman and J. Karle, submitted to Phys. Rev. 
(b) J. Karle and H. Hauptman, submitted to Acta Cryst. 
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shows that for values of a usually used, (0.005-0.010), 
the exponential function is practically a delta-function 
so that f(r) is merely a smeared or broadened ver- 
sion of D(r). The main contribution of the integrand 
of Eq. (7) is in the vicinity of p=r. 

The radial distribution curve yields three kinds of in- 
formation; the equilibrium interatomic distances, an 
evaluation of the vibrational motion between pairs of 
atoms, and the coefficients for the molecular scattering 
expression which are related to the areas under the 
peaks. The structures of many molecules can be com- 
pletely determined from radial distribution calculations 
alone. If several interatomic distances in a molecule 
fall close together, however, the information obtainable 
from a radial distribution curve is restricted and should 
be supplemented by additional data which can be ob- 
tained from computed intensity curves. The formula 
to be used is 


Ta(s)=¥ ¥ cs; exp(—3*/4hss) sinsr:;/sriz (649), (9) 


1 j=l 


where the values of ¢;;, 4;;, and r;; for various models 
are chosen to be consistent with the data available from 
the radial distribution curve. Equation (9) has been 
obtained by James’ by integrating Eq. (3), using a 
probability distribution function which is derived on 
the basis of Hooke’s law forces. 


‘PROCEDURE 


The diffraction photographs were taken with an 
electron diffraction camera which has been constructed 
recently in our laboratory. The accelerating voltage 
usually used is about 40,000 volts and is regulated to 
about one volt in 10,000. The voltage can be read con- 
tinuously on a Type K potentiometer. The voltage ob- 
tained from the measurement of diffraction rings on a 
gold photograph agrees to within 0.1 percent with that 
read on the potentiometer. The electron beam is fo- 
cused by a focusing grid and two magnetic lenses. The 
beam passes through a defining aperture of 0.1 mm di- 
ameter and is essentially parallel between the aperture 
and the photographic plate. 

A rotating sector is incorporated in the specimen 
chamber. The sector is attached to a ball-bearing race 
six inches in diameter which is rotated by a synchronous 
motor using a friction drive. The motor is in the speci- 
men chamber and its operation does not affect the 
vacuum nor the electron beam. When the sector reaches 
a speed of 1300-1400 r.p.m., the motor is turned off 
and pulled away from the race. The speed decreases 
to about 800 r.p.m. in five seconds during which time an 
exposure is made. An s* sector which is cut to operate 
between s values of one and 35 has been used so far. 
The portion between s values of zero to one acts as a 
beam stop. In order to facilitate the centering of the 
sector with respect to the electron beam, a speck of 
fluorescent material is placed in a 0.1-mm diameter 
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hole bored in the center of the sector. Photographs of 
the beam using this arrangement show that there is no 
correction for fogging required for the exposure times 
ordinarily used. 

The diffraction patterns obtained on contrast lan- 
tern slides are traced on a Leeds and Northrup micro- 
photometer. During the tracing they are rotated 
rapidly® in order to average out any irregularities 
caused by graininess. Care is taken to set the zero of the 
microphotometer on an unexposed portion of the plate. 
Traces obtained by this procedure are shown in Fig. 1. 
The microphotometer readings are calibrated to read 
electron intensities by taking several series of trans- 
mission photographs of evaporated gold on Formvar 
for various exposure times. The results of such a pro- 
cedure are easily interpreted since the reciprocity law 
holds for fast elections,’® that is, the exposure is equal to 
the product of the intensity and the time. The gold 
photographs are taken using the rotating sector and are 
rotated during the tracing in order to duplicate as 
closely as possible the procedure used with the gas 
photographs. A typical calibration curve is illustrated 
in Fig. 2. 

The microphotometer readings are converted to 
intensity readings which are plotted against the variable 
s and a smooth background is drawn through the 
oscillations.". It usually is easier to draw the back- 
ground line beyond s values of 5 or 6 if the experimental 
intensity curve is multiplied by s or s* to accentuate 
the oscillations about the background. In this respect 
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Fic. 2. Typical calibration curve for converting microphotometer 
readings to electron intensities. 


* Karle, Hoober, and Karle, J. Chem. Phys. 15, 765 (1947). 

(a) A. Becker and E. Kipphan, Ann. d. Physik (5) 10, 15 
(1931). (b) W. Bothe, Zeits. f. Physik 8, 243 (1922). 

4 Uncertainties in the position of the background line affect 
mainly the shapes but not the positions of the maxima in the 
radial distribution curve. Hence, fairly reliable equilibrium dis- 
tances in a molecule can be obtained from the first computation of 
a radial distribution curve and an intensity curve based on these 
distances may be computed. The background line may then be 
corrected if necessary to correspond to the position of the zero 
axis in the intensity curve provided that no sudden changes occur 
in curvature of this line. This procedure may be repeated until 
it is evident that no additional improvement is possible. The 
improvement in the position of the background line should ordi- 
narily be accompanied by a decrease in the magnitude of the 
negative portions of the radial distribution curve. Ideally, the 
radial distribution curve should be positive. 


ELECTRON DIFFRACTION STUDIES 














i. 4 +. 
fe) 


v 5 10 is 20 25 ~o Ss 


Fic. 3. The upper curve represents the intensity obtained from 
a fairly dark photograph of CCl, taken with an s* sector. The 


lower curve represents the intensity from the u curve multi- 
plied by s*. Background curves are drawn through the oscillations. 


our experiments would probably have been aided by 
the use of a sector cut to a power of s higher than 2. 
Figure 3 illustrates how the background line was drawn 
for an intensity curve from a photograph of CCl, whose 
useful density range covers values from seven to 27. 
The: molecular scattering is obtained by dividing the 
experimental intensity by the background curve. This 
procedure not only eliminates the background scatter- 
ing but replaces by constants, except at low s values, 
the variable coefficients in the molecular scattering ex- 
pression. It is possible to measure the intensity of 
molecular scattering to much higher accuracy than 
may be implied by the inaccuracies in determining the 
absolute intensity. The inaccuracy in absolute measure- 
ment may be greater than the magnitude of an oscilla- 
tion at the larger scattering angles. However, molecular 
intensity measurements concern the magnitude of an 
oscillation relative to a background line. The accuracy 
is therefore determined by the amplitude of an oscilla- 
tion rather than by the measurement of total photo- 
graphic intensity. Five photographs ranging from light 
to heavy exposures were used to obtain the molecular 
scattering curves. The uppermost curves in Figs. 4 and 
5 are the experimental molecular scattering curves for 
CCl, and CO;. The upper curve in Fig. 5 shows the 
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Fic. 4. The three curves represent the experimental molecular 
scattering, the calculated molecular scattering, and the radial 
distribution for CCl,. The vertical scales for the molecular scatter- 
ing curves are the same 


spread of experimental points for CO2. The magnitudes 
of the maxima and minima in this experimental curve 
are known to within five percent. The characteristic 
asymmetry in the maxima and minima is well defined. 

A radial distribution curve is calculated using the 
experimental data. The integral in Eq. (6) is replaced 
by a summation where the s interval is 0.2. For the CCl, 
and CO, computations, maximum s was 25 and @ was 
chosen to be 0.0057. The calculations were performed 
with the aid of punched cards and I.B.M. machines.” 
In the experimental intensity curve there is no useful 
data below s values of two and the region below s 
values of three or four does not represent a molecular 
scattering curve with constant coefficients. Hence, the 
res.on between s values of zero to three or four is re- 
placed by an intensity curve computed with constant 
coefficients from an assumed model. If this model is 
inconsistent with the results of the radial distribution 
calculation, the calculated portion of the intensity 
curve is adjusted and a new radial distribution curve is 
computed. With this method of successive approxima- 
tions, relatively few calculations need be made. The 
radial distribution curve should consist of a series of 


2 Shaffer, Schomaker, and Pauling, J. Chem. Phys. 14, 659 
(1946). 
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Fic. 5. The three curves represent the experimental molecular 
scattering, the calculated molecular scattering, and the radial 
distribution for CO. The dots in the upper curve are the experi- 
mental points obtained in five independent experiments. Many of 
the points lie on the curve itself. The vertical scales for the molecu- 
lar scattering curves are the same. 


positive peaks, since the scattering from the electronic 
shells of the atoms has been eliminated. The radial 
distribution curves for CCl, and CO, are illustrated in 
Figs. 4 and 5.* 

Each peak of the radial distribution curve is fitted 
with the function B;; exp[—i;(r—1;;)*/(40h:;+1)], 
Eq. (8), to obtain A,; and r;; values. Only the upper two- 
thirds of a peak is used in order to minimize the effect 
from small oscillations in the base line. Relative values 
of the coefficients c;; which are used in the calculated 
intensity curves can be obtained in two ways; (a) each 
ci; value is proportional to the area under its corre- 
sponding peak multiplied by the equilibrium distance 
associated with that peak, or (b) the maximum value 
of the peak, B;;, is proportional to ¢;;(4;;)*/r;;(4ah;;+1)*. 
Using the experimental h,;, r;;, and c;; values, theoreti- 
cal intensity curves were computed for CCl, and CO, 


* Note added in proof.—Recent investigations in our laboratory 
on more complicated molecules have demonstrated that extra- 
neous oscillations such as those which appear in the f(r) curves 
of Figs. 4 and 5 can be considerably reduced by readjusting the 
background line in the total intensity curves. It has also been 
found that a peak in the radial distribution curve formed from 
several interatomic distances may be decomposed into its com- 
ponents which yield directly the detailed structure of the molecule. 
Computed molecular intensity curves are used mainly to establish 
the range of uncertainty. 
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using Eq. (9). As expected, these were found to agree 
very well with the experimental curves (Figs. 4 and 5). 
In this instance, the calculated intensity curves merely 
illustrate the consistency of the procedure. On the 
other hand, when there are structural features that are 
not resolved in the radial distribution curve, more 
emphasis will have to be placed on theoretica! intensity 
curves. 


RESULTS 


The experimental results are listed in Table I. For 
CO:2, the C—O distance compares favorably with the 
spectroscopic result of 1.163A for the equilibrium value. 
The experimentally determined values for the O—O 
distance is not exactly double the value for the C—O 
distance. This deviation, however, falls within the. ex- 
perimental uncertainty. The average amplitude of the 
vibrations of the C—O bond was calculated from spec- 
troscopic data using the expression‘ 


P=1/2h= (mrt u2)/2a, (10) 
where a is the force constant between the C and O 
atoms and uy is the average energy of an oscillator given 
by Planck’s expression, n= hv/2+hv/(exp(hv/8T)—1). 
The expression for the O—O distance is 


P=y1/a. (11) 


The numerical values used were a=15.5X10° dynes 
cm-!, v;=4.0X10" sec, and »,=7.0X10" sec..8 
The electron diffraction value for (P)! for the C—O 
distance agrees well with the spectroscopic value. The 
difference between the electron diffraction and the 
spectroscopic values of (/?)! for the O—O distance may 
be attributed mainly to the fact that the O—O peak 
in the radial distribution curve is small which makes the 
determination less accurate. In addition, it should be 
noted that the spectroscopic value is computed from a 
simplified theory (Eq. (11)). 

There are no spectroscopic values for the interatomic 
distances in CCl, available for comparison. Our values 
may be compared with other electron diffraction de- 
terminations. Finbak and Hassel*® obtained 1.770A 
and 2.876A for the C—Cl and Cl—Cl distances using 
their sector method. Pauling and Brockway“ obtained 
1.76A and 2.87A for the C—Cl and CIl—Cl distances 
using the visual method. The spectroscopic values for 
the vibrational motion between the C—Cl and Cl—Cl 
pairs listed in Table I have been computed by James’ 
using the central force model. Rosenthal'® has shown 
that CH, is the only tetrahedral molecule to which the 
central force model is applicable. By applying a general 
force field to tetrahedral molecules, Rosenthal'® ob- 
tained a value of 3.6X10° dynes cm™ for the C—Cl 


3G. Herzberg, Infra-Red and Raman Spectra of Polyatomic 
Molecules (D. Van Nostrand Company, Inc., New York, 1945), 
. 173. 


4 L. Pauling and L. O. pane, J. Chem. Phys. 2, 867 (1934). 
30 


16 J. Rosenthal, Phys. Rev. 46, (1934). 
16 J. Rosenthal, Phys. Rev. 49, 535 (1936). 
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CC 
1.770+0.010 
2.877+0.020 


0.041+0.005 
0.054+0.005 


0.055 
0.071 


CO; 
1.162+0.010 
2.310+0.020 


0.034+0.003 
0.040+0.007 


0.034 
0.029 








force constant in CCl. This is larger than the value 
2.0X 10° dynes cm used by James and implies that the 
amplitude for the vibrational motion obtained by James 
is too large. Better agreement between the electron 
diffraction and spectroscopic evaluations of the vibra- 
tional motion can be expected therefore from a recalcu- 
lation of the spectroscopic values based on the more 
general force model for CCl. 

It should be mentioned that even very small vibra- 
tions like those found in CCl, and CO, damp the in- 
tensity of molecular scattering considerably. For CCl, 
the intensity at s equal to 25 is 60 percent less than it 
would be for a rigid structure while for CO, the in- 
tensity is 35 percent less. 

The c;; values were obtained from the radial distribu- 
tion peaks using the methods given above. For CO; the 
value found for the ratio cc_o/co-o is 1.50+0.02. This 
value may be compared to the ratio 2ZcZo/ZoZo which 
also is equal to 1.50. Zc and Zo are the atomic numbers 
of carbon and oxygen. For CCl, the experimental values 
of the ratio ¢ci_ci/¢c_c) is 4.22+0.04 whereas the ratio 
3ZeiZci/2ZcZci equals 4.25. These figures indicate 
that the scattering factors of the atoms involved have 
about the same shape and that their variability with re- 
spect to the angle variable s can be eliminated by divid- 
ing the total scattering by the background scattering. 

An approximation was made in Eq. (5) by replacing 
the variable r in the denominator by the equilibrium 
distance r;;. The effect of this approximation on the 
determination of the equilibrium distances is very 
small. The correction is equal to 1/2h;7;; and amounts 
to not more than 0.001A for the CO, and CCl, mole- 
cules. There is a negligible effect on the shape of the 
peak. 

DISCUSSION 


The results on CCl, and CO, indicate that accurate 
intensity data can be obtained from electron diffraction 
investigations of gases. This is borne out by the results 
in Table I and also by the fact that the radial distribu- 
tion curves are essentially positive. Relatively little 
error in the experimental intensity can yield large nega- 
tive fluctuations in the radial distribution curves. 
By the procedure described above the magnitude and 
shape of the various features of the diffraction pattern 
can be reliably determined, thus minimizing consider- 
ably the uncertainties which are bound to arise in the 
visual examination of the diffraction photographs. The 
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Fic. 6. The solid lines are the experimental scattering curves 
obtained with an s* sector for CC and CO:. The dotted lines 
represent the background scattering calculated from tabulated 
values and multiplied by s*. 


visual method can be used with confidence to de- 
termine the equilibrium distances in a molecule if the 
diffraction pattern consists of well-defined maxima 
and minima and if the determination can be based 
mainly on the positions of these maxima and minima. 
When the determination of a molecular structure de- 
pends upon the relative intensities of the maxima and 
minima or upon special features such as shelves or 
asymmetries, the results from the visual method may be 
considerably less reliable. 

There is sufficient difference between the procedure 
developed by the Norwegian school and the one pre- 
sented here to merit some comment. In their analysis 
of the diffraction data, the molecular scattering curve is 
obtained using the expression” 


L(Z.- FP +LSi 





» (12) 


s‘ 


where P,(s) is the density reading on a microphotometer 
trace of a diffraction photograph, P,(s) is the back- 
ground curve drawn on this trace to represent the 
atomic scattering, and [> (Z:—F,)*+>.S;]/s* is the 
coherent plus incoherent atomic scattering where Z; is 
the atomic number, F; is the atomic scattering factor 
for x-rays, and S; is the incoherent scattering function 
tabulated by Bewilogua.'’ Our procedure differs in that 
the molecular scattering function with constant coeffi- 
cients is obtained from the expression 


T,(s) 
Ia(0)=K| -1| 
I,(s) 


‘TL. Bewilogua, Physik. Zeits. 32, 740 (1931). 


(13) 
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where the photographic density readings have been 
converted to intensity readings. In addition, we do not 
reintroduce the variable scattering factors once they 
have been eliminated by dividing by the background. 
The disadvantage of using Eq. (12) for the molecular 
scattering function is that a proper Fourier inversion 
cannot be made on an intensity function whose coefh- 
cients (corresponding to the c,; in Eq. (4a)) are func- 
tions of s. The use of Eq. (12) in a Fourier transform 
affects the accuracy of determining the equilibrium 
distances in a molecule very little. However, it is not 
possible to use such a curve for the evaluation of the 
vibrational motion in a molecule. 

The background scattering merits some discussion. 
It has always been considered that the background 
scattering in electron diffraction patterns of gases can 
be represented by 


(14) 


16(s)=(K/s\(X(2:— FP +LSi]. 


The results of our experiments are shown in Fig. 6 
where the dotted lines represent Eq. (14) calculated 
from tabulated values and multiplied by s*, and the 
solid line is the experimental scattering obtained with an 
s* sector for CCl, and CO,. Theoretically, the solid line 
should oscillate about the dotted line. The same type of 
discrepancy is apparent from the examination of pub- 
lished microphotometer traces by the Norwegian school? 
and by Yearian and Barss.'* Experiments on the scatter- 
ing from free atoms have been performed and are being 
continued to establish with certainty that the shape of 
the background is not a result of some error in the 
technique of performing the diffraction experiment. 
The consistency obtained in our experiments to date in- 
dicates that the theory of the background scattering 
may be incorrect. It is obvious, however, from the 
procedure described in this paper that it is not neces- 
sary to have a theoretical expression for the background 
scattering for the determination of molecular struc- 
tures. 
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for performing electron diffraction computations. We 
also wish to thank the Naval Ordnance Laboratory for 
making available their mechanized analysis facilities 
and the Washington Office of the International Business 
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The structure of p,p’-dimethoxybenzophenone has been determined from X-ray single-crystal 
data. The space group is P2,/a and there are eight molecules in the unit cell. The cell parameters 
16-43, b 16-03, ¢ = 9-62 A, B 100° 15’. The phases of the largest 270 normalized strue- 
ture factors, Ey, were determined directly by probability methods (Hauptman & Karle). A Fourier 
map with these Ey as coefficients revealed the structure, which was improved by the least - 
squares method. A line connecting the methoxy oxygens within a molecule is almost parallel to 
the 6 axis and a second molecule is shifted by approximately 4c from the first, except that the 
corresponding benzene rings of the adjacent molecules are twisted in opposite directions. The 
methyl carbons lie almost in the planes of the adjacent benzene rings. 

Also included in this paper is a discussion of the limitations of 2, for P2,/a and the theorctical 


are a 


hasis for the use of an E map. 


1. Introduction 


By this time several structures have been determined 
using the direct probability method presented in 
Monograph I (Hauptman & Karle, 1953). Christ, 
Clark & Evans (1954) have found the structure of 
colemanite and Christ & Clark (1956) have found the 
structure of meyerhofferite. Bertaut & Blum (1956) 
have derived the forms for 2,, 2, and 2, of Mono- 
graph I appropriate to the space group Bbmm, and 
have applied them to solving the structure of Ti,CaO,. 

Vand & Pepinsky (1953) purported to demonstrate 
the failure of the methods based on the joint prob- 
ability distributions of Monograph I by means of a 
simple planar four-atom structure. However, Bertaut 
(1955), using the probability formulas appropriate to 
the correct space group, p4m, easily obtained the true 
structure, thus nullifying the claim based on this four- 
atom problem. 

The X-ray structure analysis of p,p’-dimethoxy- 
benzophenone, 


O 


\ / 
\—-C~ 


CH,-O~ } ~O-CH,, 


the most complex structure thus far studied by direct 
probability methods, is described in this paper. Signs 
for the largest normalized structure factors (Z) were 
determined directly from the magnitudes of the in- 
tensities alone by the probability methods of Mono- 
graph I. These signs were used to compute a three- 
dimensional E map, which revealed the structure 
immediately. The coordinates obtained from this map 
were then improved using a three-dimensional least- 
squares procedure. A note concerning this structure 
has been published by us (Karle, Hauptman, Karle & 
Wing, 1957). 


2. Experimental measurements 


The crystals used in this investigation were grown bv 
Mr C. E. Miller of Bell Telephone Laboratories. Dr 
S. Geller suggested that we apply the probability 
methods to this material and supplied us with crys- 
tallographic data and preliminary diffraction photo- 
graphs. 

The crystals were colorless, opaque needles with the 
c axis as the needle axis. The dimensions of the unit 
cell were obtained from rotation and Weissenberg 
photographs. The crystallographic data are: 


a = 16-43, b = 16-03, c = 9-62 A, 8 = 100° 15’, 
space group P2,/a,Z = 8, 
o(X-ray) = 1-289, (physical) = 1-259. 


There are 30 carbon atoms, 6 oxygen atoms, and 28 
hydrogen atoms in general positions in an asymmetric 
unit. 

Equi-inclination Weissenberg photographs were 
taken about the three crystallographic axes, using 
nickel-filtered copper radiation. Five exposure se- 
quences were made of each layer. The intensities 
were estimated visually, using a calibrated comparison 
strip, and corrections were made for Lorentz and 
polarization factors, obliquity, «,, «, resolution and 
spot size (Wing & Birks, 1954). 5527 independent 
reflections were measured, almost the entire copper 
sphere, of which 1261 intensities were observed to be 
zero. The data were corrected for vibrational motion 
and placed on an absolute scale by means of a K curve 
(Karle & Hauptman, 1953; see also Wilson, 1949), 
thus giving Fj. From these, the magnitudes of the 


normalized structure factors were computed by means 
of 


\ 
2 215. ~F2 
Ey, File Sf: 
j=1 
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where, for spaee group P2,/a, ¢ ~ 2 when h is AOl or 
Ok) and ¢ | otherwise, N is the number of atoms 
in the unit cell, and fj is the atomic scattering factor 
of the jth atom. 

The theoretical probability distribution of a struc- 
ture factor in a centrosymmetric crystal predicts 32% 
of all |B 1, 5%, of all |B] > 2, and 0-3% of all 
E 3. The actual distribution of |Z£| for this crystal 
was found to be 25% greater than 1, 6% greater than 
2 and 1-4% greater than 3. Theoretically, for a centro- 
symmetric crystal the following averages occur: 


(\E\) = 0-798, (B®) — 1-0 and <|E?—1)) = 0-968. 


For this crystal the values for these averages were: 


(\E\) = 0-675, (h*) = 0-986, and <«j£*—1|) = 1-17. 


3. Determination of the phases 


Sign determination by the probability formulas in 
Monograph I* requires only the knowledge of the 
space group, and the magnitudes of the normalized 
structure factors. The chemical composition of the 
unit cell, although useful, is not required for the struc- 
ture determination. To facilitate the determination of 
the phases, the Zax: were arranged in order of descend- 
ing E magnitudes in groups where the indices were 
gig, 090, ggg, ggu, guu, gug, ugu, ugg, uug, and uuu 
(g = even, u = uneven). The £*—1 values were also 
listed. It was immediately apparent from the £ values 
that the reflections with / even were considerably 
stronger than those with / odd. Furthermore, in the 
groups ggg and ugu the reflections with k=4,8,12,... 
were much stronger than those with k = 2,6, 10, ... 
and in the groups ugg and ggu all the strong reflections 
had k = 2, 6, 10, .... There was no particular pattern 
for k odd. Since an ordered distribution of normalized 
intensities (Z*) with respect to the k index existed for 
this crystal, it was necessary to examine the effect it 
would have on the phase-determining formulas. 

The square of the normalized structure factor for 
P2,/a for the case of equal atoms is 


\ 


Ei, = Fi,,/ Xf? 


16 “9 cos? cos? 


220 (har; + lz;) 3 2a ky; 
sin 


N ;-, sin? 


16 */ cos 


N *~. sin 


i+ 


cos 
2a(ha;\lz;) . 2a(hx, 4-lz;) 
sin 


cos cos 
x Qaukys —. 2cky, , (3-1) 
sin 


sin 


where the cosine terms apply to the case where h+k 
is even and the sine terms to h+ odd. Using various 


* Since the solution of this erystal structure, more powerful 
phase-determining formulas have been derived (see Haupt- 
man & Karle, 1957, 1958 and Karle & Hauptman, 1957, 1958). 
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trigonometric relations, (3-1) can be rewritten in the 
form 
] 


Eva! _ “ 1)’ yi 


ee (3-2) 


(—1)"**(B3,,—1) Buwat(—ltR, = (3-3) 


l 
Ni 


where 


4 4 

R= — SX [(-1)''* 4+ cos tar(hay + lzy)] cos 4rky; 
N jek 

4 v4 

t= ¥ {(—1)"'* cos 2afh(z;4-7,-)4 Uz 

N p+ 


+ cos 2a[h(x,—:x,-) + U(z,—z,)]} 


z,)| 


x {(—1)'"* cos 2ak(y;+ y,) + cos 2ak(y,—y,)} . 


Upon averaging over the index &, expressions (3-2) 
and (3-3) become 

(Bin—Ve = CR + (Red + Ra) (3-4) 
and 


’ l 
((—1)* (Bia), = Ni Bx out Roi + (Rs (Rey 


(3-5) 
where we find, for example, 


4 
(R,, = Y 


v4 
= cos 2a[h(x,—2,) | Uz,—2,)], 


(33-4) 
+r 

4 Va 
(Ry = 5 X 008 2afh(r, 4 


+’ 


ay) 4-U(z,-2,)], (3-7) 
and the indices » and »’ range over those values of 
j=1,2,...,N for which y,—y, = 0, (y, + 0 or 4). 
In other words, the subscripts », »’ refer only to those 
pairs of atoms whose y coordinates (+ 0 or }) are the 
same. Similar expressions are obtained for (R,), and 
(R,>, with y,+y, = 0 (y, + 0 or $) and for (R,), and 
(R,), with y, = y, = 0 or }. In the general case where 
y, + +y, and y, + 0 or }, then all the remainder terms 
<R,), reduce to zero. Consequently, in the general 
case, the value of (3-4) is zero, and (3°6) becomes the 
phase-determining formula 2, for P2,/a (Hauptman 
& Karle, 1953). Obviously, if many pairs of atoms 
in the asymmetric unit have the same y coordinate or 
are in the planes y = 0 or }, the remainder term in 
(3-5) can become significant and the signs for E,, , ., 
determined from 2, may not be reliable. 

An estimate of the feasibility for using 2, in P2,/a 
can be made from the value of the quantity 
(i Eiu— sly.» which is first averaged over k for 
fixed. h, 1 and then over all h, l. If the y coordinates 
of the atoms are not equal to 0 or } or to each other, 
then the above quantity should have a value near zero. 
Deviations from zero are expected since a finite 
number of data are used in the averages. In this 
crystal, however, the value of the average was 0-608, 
two to three times the value to be expected from the 
size of sample used if the magnitudes of no two y 
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coordinates were identical. Thus the conditions for the 
exact validity of 2, for P2,/a were not fulfilled. 
Therefore, the use of this formula would require the 
establishment of a high level of rejection, i.e. a sign 
would be accepted only if the value of 2, were un- 
usually large when compared with its standard devia- 
tion. As it turned out, the values of 2, for P2,/a were 
sufficiently large for only three of the E,,, the 
niagnitudes of which were too small to be useful in 
the early stagés of the phase-determining procedure. 
Accordingly\£, for P2,/a was not useful, and it was 
decided to treat the crystal as if it were in PI. 

A similar analysis can be made for the 2, formula 
for Ey», 9 in which case atoms having both the same 
r and z coordinates introduce a limitation. For this 
crystal the value of ¢|(Ej,,—1),,;)|), was 0-221, a value 
several times larger than that expected from the 
number of terms in this average. Accordingly, only 
one sign was determined using this formula, Ey 1 9 = 

2-41. As it happened this sign was not useful in the 
earlier stages of the phase determination, but later on 
played the role of corroborating phases obtained 
otherwise. 

The 5, formula in space group P1 is not subject to 
the same limitations as 2, in P2,/a. Therefore the 
initial phases from which the rest of the phases were 
obtained were determined using 2, for P1, 


sE., ~ 3(Ej,- l), (3-8) 


where s symbolizes ‘sign of’ and ~ means ‘probably is’. 
Here a phase is determined only from the magnitude 
of one other reflection. In order to have a high prob- 
ability that the phase indication is correct, it is neces- 
sary that both E,, and Ej,—1 have large magnitudes. 
Since the largest negative value that Z?—1 can have 
is —1l, whereas the positive values can be much 


Table | 
i. (Bey) 


(a) (b) 
O-87 
0-04 
0-90 
0-80 


0-97 
~ 1-00 
~ 1-00 

0-84 
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greater than +1, only positive signs are obtained with 
2, in a structure of this complexity. For this crystal, 
signs for four reflections were accepted with this 
formula. They are listed in Table 1. 

P,,,(E.y) indicates the probability that a particular 
sign is positive. Values in column (a) were computed 
from the hyperbolic tangent formula of Woolfson 
(1954) and those in column (4) from (3-29) in Mono- 
graph I. Both formulas give approximate values of the 
probability but they are sufficiently accurate to 
indicate whether a sign should be accepted. A more 
accurate value of the probability can be obtained by 
taking higher-order approximations in deriving the 
joint probability (see (3-03) and p. 43 in Monograph I). 

Next, formula 2, was used to expand the table of 
signs, 

(3-9) 


8B, ~ 8 S Ex (Eysu—1). 
k 


In this formula, the signs of several E,, must be 
known in order to determine the sign of an E,., 
Again, only reflections with large E values are ex- 
amined so that the probability of assigning a correct 
sign is high. An example of applying 2, to find the 
sign of E;5 , , is shown in Table 2. For this application 
it is convenient to have a numerical listing of all the 
indices in dictionary order with their |Z| and E?—1 
values, in order to look up the necessary Ej, . , —1 terms. 
From the 10,4, 4 term alone, the P,(Ej., ,) is 0-88. 
When the other Z,, with known signs were included 
in the summation, the P,(Zj5..,) rose to almost 
certainty. 

Only those signs were accepted for which the cor- 
responding probability was high. Of course, there 
were some £,,, for which 2, was indeterminate. After 
quite a few + signs had been assigned, it was possible 
to get — signs from 2, with a high probability. An 
example is illustrated in Table 3. 

Since E?— 1 cannot be greater than — 1 in the negative 
direction, many negative terms are needed to get a 
large P_(E£.,). The one term alone in the example 
gave P_(£.,) = 0-61. It was found, however, that a 
great majority of the contributing terms gave an 
indication of a negative sign. 

When about 15-25 signs had been found with J,, 


Table 2 


=1288 
= 10 44 10 


(Ey = 5-67) 


44 04 
22 165 


4 
= 166 26 
= +(-04 + 3-45 — 1-00 + 13-3: 
sEsy ~ +3-18( + 0-044 3-45— 1-004 13-35) +. 
Table 3 
(Esp 3-00) 
1288 
683 
—0-92 


1288 
605 
— 0-86 


iz & 8 
603 
0-09 


sE2y ~ +5-67(—0-79—0-86—0-92—0-99) +... 
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they were reiterated in 2, to check for internal con- 
sistency. At this point there were a sufficient number 
of known signs to use p me 

sE,~s 2 EyEy,y - (3-10) 

k 

In this formula, the signs of both FE, and 2£y,, need 
to be known in order to obtain the sign of #,. Again, 
only the terms with the largest Z values were used to 
establish a basic set of signs. The signs determined by 
+, were confirmed by using 2, whenever possible, e.g 
the signs of 4, 4, 10 and 8, 4, 2 were assigned from 2, 
and they combined to give the sign of 12, 8, 8: 


4410 
84 2 


Ey = +3-19 
Ensk = +2°35 
128 8 Ey = 5-67 


sE¥qs, a3 ~ (+3-19)(+2-35) 
The probability that Ej. ,,, is + is 0-99 and is in agree- 
ment with the positive result from 2,. As many com- 
binations of 2, as possible should be used for every 
sign determined. For the Z, with the largest magni- 
tudes, there should be extremely few, if any, in- 
consistent sign indications with 2,. As one proceeds to 
smaller magnitudes of Z,, then some contradictions 
will occur and, at times, the sign may be indeterminate. 
So far, the only signs determined have been in the 
group ggg (including 0g0 and g0g). It seems to be a 
good policy to determine the signs of only the largest 
Egqq and then to proceed to the iargest Ey in the other 
groups. Otherwise a false pattern of signs may be 
established which will be found to give large in- 
consistencies in the other groups. In order to begin in 
the other groups, it is necessary to arbitrarily assign 
three signs which fix the origin (Monograph I). 
Ordinarily, the largest EZ, are chosen for the assigned 
signs. In this crystal, the ugg group had large E values, 
hence the largest, Zi; ,, = 7-17, was made positive. 
The next largest ugg was 11, 6,6 and 2, with the 
known signs from the ggg group was used to determine 
its sign. There were two contributors: 


i364 


202 


7-17) 13 6 
2 


+3-69) and 21 


1166 6-53 Il 66 6-5! 


8H ii.6.6 ~™ (4 7-17)(+ 3-69) + (—7-17)(— 2-29) 
With the signs for 13, 6, 4 and 11, 6, 6, signs for many 
more Eugg were determined with 2,. Arbitrary signs 
were also assigned to the largest wug and ggu reflec- 
tions and others in these groups were determined with 
2,. Finally, the signs in the four remaining groups were 
dependent upon the signs already obtained. For in- 
stance, in the gug group, combinations such as ugg+ 
uug, ggu+guu, ugu+uuu, and ggg+gug were used 
in 2). 

Two hundred and seventy signs with |Z| > 2 were 
assigned using 2,, 2, and 2,. Some cross checking for 
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consistency in groups other than ggg was performed 
using 2. Since the sign determination was carried out 
by hand computing, it was not feasible to apply 2, 
very extensively. Of the 270 phases, all but one sub- 
sequently proved to be correct. 

In a set of 1732 2, combinations which were used 
to determine the signs, there were 116 contradictions 
to the correct sign. These contradictions did not 
usually interfere with the sign determination since 
there were several stronger terms indicating the correct 
sign. Table 4 lists the number of exceptions to L, anil 
the corresponding E,£,E,/N* product associated with 


Table 4 
E,E,E,/Nt Exceptions 
5-0-11-0 0 
4-0— 4-9 l 
3-0- 3-9 9 
2-0— 2-9 20 
1-0O- 1-9 77 


them, where N is the number of atoms in the unit ceil. 
The exceptions are related to E,E,E,/N?* since this 
quantity occurs in the probability expression for 2, 
and since it is in a normalized form applicable to any 
crystal. 

The largest exception had a value of 4-8 and it was 
the only exception out of 144 terms in the range 
4-0-11-0. The exceptions followed an expected distri- 
bution in that the number was relatively small for 
large E,E,E,/N* values and larger for 
E,E,E,/N* values. 


smaller 


4. Structure 


The 270 reflections (6% of the non-zero data) for 
which signs had been determined were used to com- 
pute a three-dimensional EZ map. This is a Fourier in 
which the Ey rather than the F,, are used as coeffi- 
cients. The Z map was characterized by 37 main 
peaks, two of which were fairly weak. One of the 
weak peaks was spurious whereas the other 36 peaks 
were identified with the atoms of the two molecules 
in the asymmetric unit. The average weight of the 
six O atoms was 310, the average weight of the thirty 
C atoms was 266, and the rapidly oscillating back- 
ground rarely reached +100. 

The coordinates of the atoms as read from the 
E map were subjected to a preliminary least-square 
refinement on the IBM 704, using the program pre- 
pared by D.Sayre. About 1100  structure-factor 
magnitudes sampled throughout the copper sphere of 
data (every fifth reflection ordered on sin #) were 
used. One temperature factor per atom was used and 
the 28 hydrogen atoms in the asymmetric unit were 
omitted. At present, the R factor for the 1100 reflee- 
tions used in the least-squares procedure, including all 
the zero reflections, is 22%*. The R factor for all 5527 


* At this point it became apparent that the least-squares 


procedure was diverging owing to approximations in the 
program, and the computation wus terminated. 





1.L. KARLE. H. 


Table 5. 


Molecule 


Fractional atomic parameters 


Molecule 








\tom ala y/b zie Atom ria 


y/b zie 
0 0-432 
0 225 
0 0-435 
C(CH,) 0-429 
((CH,) 0-387 
(CO) 0-292 
0-398 
0-361 
0-327 
0-329 
0-366 
0-404 
0-327 
0-278 
0-311 
0-394 
0-445 
0-410 


0-680 
0-370 
0-064 
0-751 
0-985 
0-372 
0-607 
0-595 
0-515 
0-451 
0-467 
0-545 
0-289 
0-219 
0-140 
0-137 
0-204 
0-282 


0-077 0 
0-220 Oo 0-227 
0-067 0 0-424 
0-185 C(CH,) 0-381 
0-062 C(CH,) 0-419 
0-174 C(CO) 0-294 
0-112 Cc 0-396 
0-225 ( 0-311 
0-239 ( 0-282 
0-149 ( 0-328 
0-040 ( 0-412 
0-013 ( 0-446 
0-135 ( 0-330 
0-116 Cc 0-327 

Cc 

( 

( 

Cc 


0-434 0-686 
0-376 
0-069 
0-763 
0-995 
0-378 
0-612 
0-608 
0-526 
0-457 
0-466 
0-544 
0-296 
0-228 
0-153 
0-146 
0-210 
0-286 


0-585 
0-731 
0-543 
O-578 
0-629 
0-678 
0-612 
0-611 
0-631 
0-651 
0-650 
0-625 
0-649 
0-730 
0-708 
0-594 
0-498 
0-531 


0-091 0-357 
0-099 0-391 
0-107 0-397 
0-137 0-367 
reflections is 26%, and 22% if the reflections read to 
be zero are omitted. It would be worth while to refine 
the structure further both with three-dimensional 
Fouriers and with least squares, using anisotropic 
temperature factors and the complete set of data. 
The present atomic coordinates are listed in Table 5. 
The two molecules in an asymmetric unit are illus- 
trated in Fig. 1, a projection along the a axis, }<2< }. 


HAUPTMAN, J. 


KARKLE AND A. B. WING 














T 


b 
t 


Le C 


Fig. 1. Projection of an asymmetric unit along the « axis 
*<z7<}. 


The three oxygen atoms in each molecule lie near the 














Fig. 2. 


Electron-density projection on (001). The pairs of superimposed molecules are separated by approximately jc. 


The coordinates of the line drawings were taken from the preliminary three-dimensional least-squares refinement. 
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202 planes. The molecules are related to each other 
approximately by a shift of z+} and a reflection 
* The carbon atoms in the methoxy 
sroups are very nearly in the same plane as the 
xdjacent benzene rings. The’ top ring in molecule I 
and the bottom ring in molecule IT are rotated about 
38° out of the planes formed by the three oxygen atoms 
in each molecule, whereas the bottom ring in molecule | 
and the top ring in molecule If are rotated about 24° 
out of the plane of oxygen atoms. 

Fig. 2 illustrates the contents of the unit cell as 
viewed down the ¢ axis. The pairs of superimposed 
molecules are separated by approximately jc. The 
coordinates of the line drawings were taken from the 
least-squares cefinement. The signs for the electron- 
density projection were derived from these coordinates 
where R = 20% for the z projection, including all 
terms. The same reflections and signs were used to 
compute the E map (Epxso instead of Frxo for the 
coefficients in the Fourier series) which is illustrated 
in Fig. 3. A comparison of the E map with the F map 
shows the much greater resolution of atoms in the 
superposed benzene rings and methyl groups in the 
E map. The contours near the centers of the atoms 
were not drawn owing to the fairly rough grid used in 
computing the series. The contours in the EZ map are 
on an arbitrary scale spaced by a value of 1, with the 
lowest contour at 2. The background oscillates very 
rapidly but does not exceed +2. 

The shortest intermolecular distances are between 
an oxygen and oxygen in methoxy groups at 3-41 A 
and an oxygen and carbon at 3-45 A. The average 
bonded distances with their average deviations at the 
present stage of refinement are listed in Table 6. 


Table 6 


Average 
distance (A) 


through y 


Average 
deviation (A) 
0-030 
0-030 
0-025 
0-028 
0-010 


Bond type 


(-C, ring -378 

(.C, non-ring 475 

(-O, aromatic “405 

C—O, aliphatic -500 

c=0 -280 
The C-O-C angles are 115° and the C,,-C-C,, 
angle in the ketone group is 123°, which is the same 
as that found in p,p’-diiodobenzophenone (Manthey, 
Plieth & Singewald, 1952). 


5. Theoretical basis for E maps 


The value of removing the effect of the electron distri- 
bution and vibration in sharpening peaks in a Fourier 
series is well known. To our knowledge this sharpening 
has been applied mainly to the calculation of a Patter- 
son series. It has been seen in this paper that it is 
worth while to compute sharpened Fourier series 
representing atomic positions using on the one hand 
only the largest HZ values as coefficients and on the 
other hand all EZ values contained within a particular 
sphere in reciprocal space. A Fourier series computed 
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Fig. 3. An E map projection on (001). The same reflections 
and signs were used to compute the E map as those used in 
the F map, Fig. 2. The increased resolution over that in 
Fig. 2 is of interest. 

from a restricted number of large Z’s is well adapted 

for use with direct methods of phase determination 

since the phases of the largest Z’s are the ones which 
are most easily determined. We present here the 
theoretical significance of such a calculation. 

We define the structure factor 


Fy, os Xn, +p, (5-1) 


where 


(5-2) 


(9-3) 


Yn, = ~ Sin tity 
jan 


N is the number of atoms in the unit cell, x is the 
symmetry number, fi, is the atomic scattering factor, 
and the values of the real and imaginary parts Xy, 
and Y,, of the structure factor Fy, are known to be 
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1, and By, respectively. The € and 9 are the tri- 
gonometric functions of the atomic coordinates listed 
in International Tables (1952) for the various space 
yroups. It was shown (Karle & Hauptman, 1954) that 
the most probable values of the coordinates xz, yx, Ze 
of those atoms characterized by the atom scattering 
factor fx coincide (to good approximation) with the 
principal maxima of 


|< —(4n, Som, Fim, - (Bu, Sit. Mary)” | (5-4) 
Peer _ y 
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my, — m! =\ \ \ Fn, day dyjtlzy 


ol ol ol 


Mm, = mM, = \ \ \ Nin, 1x; dys dz; . 
0-000 
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In the cases that either m$ or mj equals zero, m, equals 
the non-zero moment. If we take the logarithm of (5-4) 
and expand the squares, we obtain 


. I , a4 
= ve (—A},,—Arn,. + 2fin, (An, Fin, + Bu, Nin,) 
“ . £2 > : + en 
amy & Simp —JSin, (Ein, +Nin,)) - (5-6) 
As pointed out by Bullough & Cruickshank (1954), 
the main term in (5-6) is 


wm 
= Sm 
== Nja 


pool vw 72 
Me on Sim 
j=2 


(An, Ein, + Bn, 7m, (57) 


the first two terms of (5-6) being a numerical constant 
and the last is space-group, though not structure, 
dependent, and is also not expected to have an im- 
portant effect upon the maxima. The last term is seen 
to be a numerical constant for space group P1. 

Expression (5-7) is a Fourier series and we are 
interested in seeing how the coefficients of the real 
and imaginary parts are related to the normalized 
structure factor Zy,, defined by 


Fy, . 
Nia = oe 
mi) S fin,)4 — ((m3+ms) =f in)? 

jul ’ 


E X hy, * iY,, 
hy, \ 


((mS + 


If we ignore numerical constants, we have for the 
coefficients of the first and second terms in (5-7), 
respectively 


(5-9) 


-J. KARLE AND A. B. WING 


263 


Samy 


Vin 


: 9 ‘ + 4 al 
(Sim) | ES) 


The first term in the product (5-9) is the same as the 
first term of (5-10). It is essentially a constant and if 
all atoms are assumed to have the same shape, it is 
an cxact constant. Now consider the second terms in 
(5-9) and (5-10). The numerators are the known values 
of the real and imaginary parts of the structure factor 
and the sum in the denominators differs from the sum 
in (5-8) by a single term. We may conclude that the 
coefficients in (5-7) are proportional to Z,,, and we are 
thus led to the significance of an E map, namely: The 
maxima of an £ map answer the following question 
“What are the most probable values of the atomic 
cvordinates given the values of any set of structure 
factors ?’. 


B 
by (5-10) 


We wish to thank many peeple for their invaluable 
assistance. Mr C. E. Miller of Bell Telephone Labora- 
tories prepared the crystals; Dr S. Geller of the same 
laboratory obtained preliminary crystallographic data 
and suggested that the probability method be applied 
to this crystal; Mr P. O’Hara of the National Bureau 
of Standards, Dr H. T. Evans of the Geological Sur- 
vey, and Mr R. A. Sparks of the University of Cali- 
fornia in Los Angeles performed the various calcula- 
tions for the K curve, the EZ maps and the electron- 
density projection; and Mr A. Hatch of the Service 
Bureau Corporation applied the least-squares program 
developed by Dr D. Sayre. 
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Profiles 
in Science 


The whole of science is nothing more than 
a refinement of everyday thinking. 


Albert Einstein (1936) 


Theodor Wieland 


Theodor Wieland is Professor Emeritus and formerly 
Director of the Chemistry Department of the Max-Planck-In- 
stitute, Heidelberg. He retired in 1981. Professor Wieland’s 
PhD thesis on toad poisons was done under the tutelage of 
his father, Professor Heinrich Wieland of the Chemical Insti- 
tute, the University of Munich. While at the University of 
Heidelberg in 1937, Wieland pursued vitamin research and 
succeeded in isolating pantothenic acid from fish liver; he 
studied their biosynthesis and the separation of animo acids 
by chromatography and paper electrophoresis. After World 
War II, he was an associate professor at the University of 
Mainz working mainly on peptide syntheses (mixed anhy- 
dride method, active amino acid thiol esters) and poisons 
from Amanita mushrooms. In 1951, he became Director of 
the Institute of Organic Chemistry at the University of Frank- 
furt and completed the structural and analytical chemistry 
of the amatoxins and the phallotoxins and isolated the anti- 
toxin cyclodecapeptide antamanide. It was this compound 
which brought him to Isabella Karle in 1973 when she agreed 
to investigate its structure. 
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Dr. Bernhard Witkop 


Dr. Bernhard Witkop retired from the National Institute 
of Health in 1993 after a distinguished tenure of 42 years. 
Since 1956 he was chief of the section on metabolites in the 
Laboratory of Chemistry of the National Institute of Arthritis 
and Metabolic Diseases. He is well known for developing the 
cyanogen bromide reaction, a method for cleaving proteins 
precisely in order to facilitate the sequencing of a protein. 
The technique was a key development in the current genetic 
engineering revolution since protein sequencing is needed 
before a gene can be properly cloned. Witkop and his cowork- 
ers have helped define the metabolic pathways of many 
biologically important molecuies, including serotonin, a 
neurotransmitter, and its precursor tryptophan, which is 
important for normal growth. For more than 30 years Bern- 
hard Witkop and Isabella Karle have collaborated on research 
projects. Most noteworthy are: (1) the structure of the venom 
of the Colombian arrow poison frog, batrachotoxin(in), at the 
time probably the record for establishing a complete struc- 
ture with an incredibly small crystal; and (2) the photolysis 
of n-chloroacetyl-o-methyl-l-tyrosine to an azulene is prob- 
ably the first application of X-ray diffraction to such a com- 
plicated molecule of completely unknown chemical 
structure. 
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Dr. Arnold Brossi 


Dr. Arnold Brossi is a visiting research professor of 
chemistry at Georgetown University, Washington, D.C. From 
1976 to 1991, he was Chief of the Section of Natural Product 
Chemistry at the National Institute of Diabetes and Digestive 
Kedney Diseases, at the National Institute of Health. His 
major research interests include study and synthesis of bio- 
logically active natural products and drugs useful in malaria 
chemotherapy. He has specialized in studies of antipodal 
isomers. His work, spanning five decades has had a great 
impact on medical research and has led to the discovery of 
several drugs (Tetrabenazine, Dehydroemetine, Dimecapto- 
succinic Acid, Astiban, Versidyne, Arteether and Phenser- 
ines). Besides using routine spectroscopic data for 
characterizing compounds, Brossi had early access to X-ray 
crystallography at the Naval Research Laboratory. Working 
with Isabella Karle and J.L. Flippen-Anderson, he used X-ray 
analysis to determine the absolute configuration of optically 
active compounds for elucidating the three-dimensional 
configuration of complex molecules and conformations of 
bioactive molecules. 


Two/1994 25 





Reprinted from Acta Crystallographica, Vol. 17, Part 7, July 1964 





PRINTED IN DENMARE 


Acta Cryst. B35 


1964). 17, 


An Application of the Symbolic Addition Method to 
the Structure of L-Arginine Dihydrate 


By IsaBeLia L. KARLE AND J. KARLE 


U.S. Naval Research Laboratory, Washington 25, D.C., U.S.A. 


(Received 25 July 1963) 


The crystal structure of L-arginine dihydrate was determined directly by means of the symbol 


using the complete three-dimensional data obtained from 
P2,2,2, and the 


uldition phase dete 
Cu Ihe 


rmination proce cdlure 


raciation. pace group is cell dimensions are: 


a=5-68, b=11-87, unl « 15-74 A. 


The argmine molecule is a zwitterion with the guanidyl group, 


rather than the amine group, 
'wo planes characterize the arginine molecule, one through the acid 
yroup, and the other through the extended side chain which contains the guanidyl group. The 
dihedral angle between these planes is 74 
infinite 


accepting an extra proton. 


. The arginine molecules and the water molecules each maks 


chains perpendicular to each other by hydrogen bonding. All 18 hydrogen atoms have been 


located. 


Introduction 
(H2N)oCNH(CHe)sCH(NHe2)COO-, is one 


plication of the symbolic addition method for phase 


s Avalos determination (Karle & Karle, to be published) to a 


of the amino acids essential for animal growth. It is 
also the most basic of the amino acids since, in addition 
to the «-amino group, it contains a terminal 
guanidyl group. This investigation afforded a good 
opportunity to study a biologically important amino 
acid structure and its associated hydrogen bonding. 
It also provided a valuable experience in the first ap- 


also 
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non-centrosymmetric crystal. 


Experimental 
Purified L-arginine was obtained from the Central 
Research Laboratories of General Mills, Inc. Re- 
crystallization from water at room temperature yielded 
crystals of the dihydrate which were colorless, transpa- 
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rent prisms. They were elongated along the a axis. 
Crystals kept in air lost their sharp edges and appeared 
to deliquesce, whereas those kept in a desiccator 
became opaque and cracked. The crystals kept in the 
X-ray beam however were fairly well preserved. 
Eventually they showed signs of disintegration in the 
region where the X-ray beam impinged upon them. 

Complete intensity data for Cu Ax radiation were 
collected along the a,b and ¢ axes by the multiple- 
film equi-inclination Weissenberg technique. Preces- 
sion photographs for the Ok/, hOl and hkO layers were alsu 
taken. No correction for absorption was made. 

The space group is P2,2;2; and the unit cell para- 
meters are: 


a=5-68 + 0-01, 


The computed density is 1-314g.cm-% and the 
measured density by flotation in mixed solvents was 
found to be 1-320 g.cm-. 


b=11-87 + 0-02, c= 15-74+0-02 A. 


Phase determination 


The normalized structure factor E, is defined by 


(1) 


where ¢e=2 when h is 00, O40 or OOL, ¢ =1 otherwise 
for space group P2;2,2;, and f; is the atomic scattering 
factor for the jth atom in the unit cell containing NV 
atoms. Phases corresponding to 400 of the largest |E, 
were determined by means of the symbolic addition 
procedure for non-centrosymmetric space groups. 
The procedure is an extension of the symbolic addition 
method applied to centrosymmetric crystals (Karle & 
Karle, 1963; Karle, Britts & Gum, 1964). 

The formula for starting the phase determination is 
(Karle & Karle, to be published), 


fn ~ GutPn—ee, » (2) 


where the g, are the phases associated with the 
E,| and k, implies that the average is taken only over 
those k corresponding to the larger |EZ| values. 
Formula (2) is implemented by specifying appropriate 
phases (Hauptman & Karle, 1956) in order to fix the 
origin, and by assigning additional symbols to re- 
present the phases associated with certain appro- 
priately chosen reflections. 

In space group P2,2;2;, it is convenient to fix the 
origin by assigning the phases of two-dimensional data. 
The phases were assigned as follows: 


h Th En 


3,0,10 0 3:46 
3,3, 0 —a2/2 2-17 
3,0, 1 +2/2 


o.77 


=~ dé 


With these assignments, it is immediately possible to 


obtain contributors to (2) for particular g,. For 
example, the largest guu is |E£¢3:|=2-06. We note that 


SYMBOLIKK 


ADDITION METHOD 


if k=301 and h—k=330, h=631. Since gss90= —2/2 
and 901 = +2/2, this pair contributes the value zero 
toward the average in (2) for ge. Similarly, since 
Yi0,;0=2% and the value zero may be used for ges 
at this stage, this pair contributes the value of 2 to 
the average in (2) for 93.311. |£3,3,1:| =2-05, the second 
largest uuu. Two pairs of contributors are now available 
for gq 0310( Eo310 = 1-85) as follows: 


h Th h 
3.0, 1 : 3,3, 0 —a2a/2 
3.3.11 3.0.10 % 


03.10 0.3.10 +a2/2 


Further contributors to new phases and also to ¢¢a1, 
@331u and go310 were developed by continuing the 
process of forming pairs as required by (2). In order 
to facilitate this, it was necessary to assign symbols 
to other phases. The assignments were, 


Yh En 
3-21 
2-31 
” 56 
2-31 


These reflections were chosen because they entered 
into many of the combinations for the average in (2) 
and they were associated with relatively large |Z 
values.* The symbols p, s,m must be either 0 or 2, 
but @ can be anywhere in the range —a<a<z. In 
forming the contributors to (2), maximum use was 
made of the three-dimensional reflections (hkl#0). 
Too much reliance on two-dimensional reflections may 
lead to misleading phase indications. 

In the course of application of the procedure, it soon 
became apparent that m=z. It was also apparent that 
a was near 0 or z. It was chosen to be zero, which, in 
effect, specified the enantiomorph. The specification 
of the enantiomorph for space group P2;2,2; is dis- 
cussed in detail by Karle & Hauptman (1956). 

By means of hand computation, 137 phases were 
determined in terms of 0,2, +2/2, p, and s. Relatively 
few inconsistencies were found for pairs involving 
three-dimensional reflections or a three-dimensional 
and a two-dimensional reflection. However, there were 
a number of inconsistencies for pairs composed solely 
of two-dimensional data. 

Since p and s have to be 0 or x, most phases in this 
first approximation were on the cardinal points, 
0,2 or +2/2. Some, however, were found to have 
intermediate values. For example, an approximately 
equal number of contributors to (2) might consist of 
the types -2/2+s and p+s. The average, 
4[{(—2/2+8)+(p+s)], would then have the possible 
values, 2/4, —2/4, 32/4 or —32/4, depending upon the 
values of p and s. 


* A procedure based on the multiplication of symbols, 


rather than addition, can be used if the symbols represent 
exp (ip) rather than @¢. 
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. . 4 
Fig. 1 
E, as coefficients. The 
data with |E min = 1-40, 


ire extraneous. The 


hree-dimensional E-maps for arginine 


phases were 


The following four combinations of values for p and s, 


were substituted into the set of 137 phases defined in 
terms of these symbols. Each of the four sets of phases 
so obtained was used as a basic set for the application 
of the tangent formula (Karle & Hauptman, 1956) 
> |Ey Ey «| Sin (Py +Pn—x) 
k, 


?, * =; , (3) 
> | Ey Ey__| ©O8 (ye +Pn—n) 
k, 


tan 


where k, ranges over 137 values of k. This formula is 
used to reiterate the determination of values for the 
initial set of phases and to obtain additional ones. 
Since the numerator of (3) is proportional to |Z, sing, 
and the denominator of (3) is proportional to |£),| cos@p, 
the sum of the squares would be proportional to |Z,,|?. 
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, viewed along the a axis, which were 
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(2) 


H,0(2) 
Cc —_— 


omputed with normalized structure factors, 


obtained from the symbolic addition procedure. (1) 137 data with |E) min= 1-50, (2) 200 
(3) 300 data with |E min= 1-15, and (4) 400 data with 
ontours are at equal intervals on an arbitrary scale. 


E\ min = 1-00. The peaks marked with crosses 


If this formula is calculated for many h values, the 
proportionality constant may be obtained from a 
comparison with the experimental values for |£,|*,.. 
over the same range of h. The adjustment is made 
by setting the average of the calculated |2,)*,), 
to the average of the observed |£,,|*,,.. In the course of 


equal 


the calculation those reflections were eliminated 
whose |£,|..;. were very small. An arbitrary level of 
0-3 was chosen. Reflections whose phases changed 
greatly from reiteration cycle to reiteration cycle 
were also eliminated. Usually, the tangent formula (3) 
converges in one or two iterations if no new data are 
introduced. 

The correct answer was found to correspond to the 
combination p=s=z. The other three combinations 
did not produce Z-maps which were crystallographi- 
cally meaningful. E-maps are Fourier maps employing 
the E, as coefficients (Karle, Hauptman, Karle & 
Wing, 1958). Fig. 1 shows E-maps calculated for (1) 
137 data (minimum |£|=1-50), (2) 200 data (mini- 
mum |E£|=1-40), (3) 300 data (minimum |Z} =1-15), 
(4) 400 data (minimum |Z|=1-00). It is seen that 
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calculation (3) of Fig. 1] reveals the molecular structure 
quite well. This map includes about one fifth of the 
total data collected. The peaks marked with crosses 
are extraneous. In calculation (4) of Fig. 1, there were 
no extraneous peaks as strong as the fourteen peaks 
which corresponded to the structure. 

[It is of interest to note that phase-determining 
formulas obtained in the past, which define phases 
in terms of structure factor magnitudes, could have 
been used to determine the correct values of p and x 
directly. One type of formula which could have been 
used is the 2, formula for space group P2)2)2; (Karle 
& Hauptman, 1956). The two forms of interest are 


SEensro ~ SEY (—1V)"(| Bnws|? —1), (4) 


SEo2ien ~ SS (—1)"**(|Bags!2—1), (5) 
where S means ‘sign of’. Table 1 shows the results of 
the application of (4) and (5). There were nineteen 
terms in the sum for (4) and seven in the sum for (5). 
The probability that the sign be positive for the E., 
in question, P,(E.,), is given by 


> 


P., (Ex) = +4 tanh o3! F.,|21/202*”? (6) 


where 3; is the sum contained in (4) or (5), and 


On= SZ", Z; being the atomic number of the jth atom 
1 

in a unit cell containing N Atoms. It is seen from 
Table 1 that three determinations indicate that p has 
a value of x with probabilities 0-86, 0-72, and 9-73, 
and s has a value of z with probabilities 0-85 and 0-72. 
The over-all probability from the three determinations 
that p is equal to z is about 0-97 and that for s from 
the two determinations is about 0-94, assuming that 
the determinations are independent. 


Table 1. An application of the 2) formula 
b ad P 
—7" 


0-14 
0-15 


ia 


S\ mbol ,(Eon) 


0 


Another type of phase-determining formula which 
could have been used for determining the value of p 
and s is B3,o (Karle & Hauptman, 1958). With this 
formula there were several indications that the values 
of p and s were x. The value of m was also found to 
be z. It is interesting to note that when the data were 
extended by use of a Patterson map, improved by the 
application of a positivity criterion (Karle & Haupt- 
man, 1964), the results of formula B3,.) were more 
definitive. 

In the initial set of 137 phases obtained from formula 
(2), the phases of the three-dimensiona] reflections had 
errors which were distributed as indicated by the 
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tial phases from eg. (2) 


phases refined by tangent 


formula (3) 


160 





Fig. 2. The maynitude of the difference hetween the initial 


values of phases for three-dimensional reflections and the 


final values calculated from the refined structure. 


dashed curve in Fig. 2. The symbol q@error is defined 
as the magnitude of the difference between a phase 
calculated for the refined structure and one calculated 
from (2) or (3). It is seen from the solid curve that 
application of (3) to the initial phases improved their 
values considerably. The average error for the three- 
dimensional phases decreased from 43° to 27 
use of the tangent formula. 


by the 


Refinement 


The preliminary coordinates of the fourteen atoms, 
other than hydrogen, which were obtained from the 
fourth E-map in Fig. 1, were subjected to a three- 
dimensional least-squares refinement (Busing, Martin 
& Levy, 1962) using all the 1406 data. After two cycles, 
employing isotropic temperature factors, the R index 
was 14-8%. Next, two cycles with anisotropic tempera- 
ture factors reduced the R index to 11-6%. All the 
terms were evenly weighted with the value of one. At 
this point, a difference Fourier synthesis was computed 
in order to find the hydrogen atoms. All 18 hydrogen 
atoms were located. One more refinement was per- 
formed which included all the 32 atoms. Then another 
difference map was computed, which used the 4F from 
the heavier atoms alone while the phase angles 
included the calculated contributions from the 
hydrogen atoms as well. The final R index was 10-3%.* 

The difference map showing the location of the 18 
hydrogen atoms is illustrated in Fig. 3. It is a composite 
drawing projected along the a axis. Although the 
coordinates of the hydrogen atoms have not been 
extensively refined, there is no doubt about the ap- 
proximate location of each hydrogen and the configura- 
tion about the heavy atom to which it is bonded. A 
final electren density map computed with F values 


* \ microfilm containing the observed F and the computed 


F, A and B has been deposited with the Library of Congress, 
Washington, D.C., U.S.A. The reel number is 8496. 
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Table 3. The approximate coordinates for the 
hydrogen atoms as determined from the second 
difference map 


Atom Bonded to z 1 z 


H(t) C(4) 0-908 . 0-350 
H(2) N(5) 0-847 . 0-328 
H(3) N(5) 0-938 ; 0-362 
H(4) C(6) 0-076 0-475 
H(5) C(6) 0-933 -525 0-495 
H(6) C(7) 0-366 0-450 
H(7) C(7) 0-208 -682 0-502 
H(8) C(8) 0- . )-589 
H(9) C(8) 0-258 . 0-623 
H(10) N(9) 0-52: . 0-617 
H(11) N(11) 0-72: 0-660 
H(12) N(11) “86: “475 0-705 
>) +18) H(13) N(12) 933 0-719 
=A H(14) N(12) -822 7 0-700 
H(16) H(15) W(1) 53% 0-629 
: : -— 
Fig. 3. A projection along the a axis of sections of a difference Hi a@) ws) , ° — 
map for the location of hydrogen atoms. Contours are at H(17) W(2) o-412 
intervals of 0-05 e.A~3, starting with the 0-15 e.A-3 contour. H(18) W(2) 7 0-463 


is shown in Fig. 4. This composite drawing is viewed 
along the 6 axis. 

The final coordinates and anisotropic temperature 
factors for the heavier atoms are listed in Table 2. 
The approximate coordinates for the hydrogen atoms 
are given in Table 3. 


Geometry of the molecule 


The arginine molecule is characterized by two planes, 
one through the carboxy! group and the other through 
the side chain which terminates with the guanidyl 
group. The molecule is a zwitterion with the guanidyl 
group, rather than the amino group, accepting a proton 
from the acid group. The bond lengths and angles 
are listed in Table 4 and shown in Fig. 5. 
; ’ The dimensions of the amino acid group are quite 
electron density map. Only one arginine molecule is shown. . 6 ; ~ 
Contours are at intervals of 1 e.A-3. starting with the 1-5 ‘Similar to those reported for glycine (Marsh, 1958) 
43 contour and L-lysine (Wright & Marsh, 1962). The equation 


Fig. 4. A projection along the 6 axis of sections of the final 


Table 2. The final parameters and their standard deviations 


thermal parameters are of the form T=exp [—(f,,h*? + Bogk* + Bggl* + 28 phk + 28, ghl + 2B,,kl)]. For the oxygen atoms of 
water molecules, the alternate symbols W(1) and W(2) are used in the succeeding Tables. The standard deviations are 
multiplied by 104 
Oz / z 0: 10*,, 10*B., 1088,, 1048,, 1088,, 10*f,, 

O(1) 0-3022 (11) 0-4036 0-3231 (4) 244 44 42 23 68 
O(2) 3340 11) 0-5832 (4) 0-2830 (4) 219 45 31 j 53 
C(3) 2353 14) 0-5047 (6) 0-3220 (4) 162 41 17 16 
(4) 0072 11) 0-5293 (5) 0-3724 (4) 95 35 14 6 
N(5) 0839 (11) 0-6432 6) 0-3524 (4) 136 44 23 0 

0622 13) 0-5186 (7) 0-4680 (4) 149 12 l 


+ 


| 
Sa me oe ND 8 OO tO 


( 

C7 2546 (17) 0-5995 (7) 0-5002 288 5¢ 19 
('(8) 3822 (14) 0-5478 (7) 0-5753 148 52 22 
N(9) 5532 (11) (6274 (5) 0-6118 (3) 27 14 
C(10) 7315 (12) 0-5955 (6) 0-6633 (4) 110 j 10 
N(11) 7796 (12) 0-4862 (5) 0-6760 (4) 161 36 16 
N(12) 8579 (13) 0-6739 (6) 0-7016 (4) 191 24 


— 


H,O(1), W(1)  0-5506 (13) 0-2790 (6) 0-6214 (5) 252 f 66 
H,O(2), W(2) 0-6222 (17) 0-2871 (12) 0-4356 323 50 
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Fig. 5. Bond distances and angles. 


for the least-squares the atoms 


O(1)O(2)C(3)C(4) is 


plane through 


3°1187x + 2-3040y + 12-795z = 6-0093 (7) 
and the deviations of the individual atoms from the 
plane are: 


0(1)—0-003, O(2) — 0-003, C(3) +0-008, and 
C(4)—0-002 A. 


Table 4. Bond lengths and angles 
Bond 


(3)-O(1) 
(3)-O(2) 


Angle 


-259 A 
-249 


O(1)-C(3)-O(2) 
O(1)-C(3)-C(4) 
O(2)—C(3)-C(4) 


125-6 
115-2 
119-1 
C(3)—C(4)—N(5) 110-5 
C(6)—C(4)—N(5) 110-7 
C(8)—N(9)-C(10) 123-2 
N(9)—C(10)—-N(11) 121-0 
N(9)—C(10)—N(12) 118-9 
N(11)-C(10)-N(12 120-2 


C(3)-C(4)-C(6) 108-4 
C(4)—C(6)-C(7) 114-4 
l 
l 


1(4)-—N(5) 480 
‘(8)—-N(9) 471 
(10)—N(9) -351 
(10-N(11) -340 
(10)-N(12) 322 
(3)—C(4) -547 
(4)-C(6) -542 
‘(6)—C(7) -540 
(7)-C(8) 517 


C(6)—C(7)-C(8) 110- 
C(7)-C(8)—N(9) 111: 

The standard deviation for the bond lengths is 0-008 
0-012 A and for the angles 0-9°. 


The nitrogen atom N(5) lies 0-280 A out of the plane; 
the values found in glycine and lysine were near 
0-44 A. Although the amino group does not accept an 
extra proton to form C-NHj as was the case, for 
example, in glycine and lysine, the nitrogen in the 
amino group does assume a tetrahedral configuration 


OF THE SYMBOLIC 


ADDITION METHOD 


with C(4), H(2) and H(3), and the fourth position on 
N(5) is occupied by a hydrogen bond from N(9’). 

The guanidy] group is attached to the carbon chain 
by a C-N bond of 1-471 A, whereas the three C-N 
bonds in the guanidy] group are nearly equal and 
average 1-338 A. The C-N bonds in urea were found 
to be 1-336 A (see, e.g. Sklar, Senko, & Post, 1961), 
1-34 A in nitroguanidine (Bryden, Burkardt, Hughes 
& Donohue, 1956), and 1-34 A in guanidinium bromate 
(Drenth, Drenth, Vos & Wiebenga, 1953). The N-C-N 
angles are nearly 120° and the guanidy! group is planar. 
The equation of the least-squares plane through the 
atoms N(9)C(10)N(11)N(12) is 


3°57642 + 0-1756y — 12-2262 = —5-3907 (8) 


and the deviations of the individual atoms from this 
plane are: 
N(9) —0-0005, C(L0)+0-0015, N(11)—0-0005 and 


N(12)—0-0005 A . 


The five hydrogen atoms, H(J0) to H(14), of the 
guanidyl group, — NHC(NH.,);, are also very nearly 
in the plane of the NCNN, thus making the whole 
guanidy! group planar. 

The average of the four C—C bond lengths in this 
molecule is 1-537 A and the average of the four 
C-C-C angles is 111-0°. However, the angle 
C(4)—C(6)—C(7) at 114-4° is somewhat larger than the 
other three. A similar widening of the C—C-C angle 
adjacent to the amino group has been observed in 
other amino acids and peptides. The side chain 
including the guanidyl group is fully extended and 
lies approximately in a plane. The equation for the 
least-squares plane through these atoms is 


3-66842 — 0-5873y — 11-9922 = — 5-6305 (9) 


and the dihedral angle between this plane and the one 
through the acid group is 74°. 


Packing of the molecules 


The packing of the arginine molecules is characterized 
by a three-dimensional network of hydrogen bonds. 
A drawing of the structure viewed along the a axis is 
shown in Fig. 6. The arginine molecules lie extended 
parallel to the c axis and they form infinite chains by 
means of the hydrogen bonding between two of the 
NH of the guanidy] group and the two O of the carboxy! 
group. The guanidy] and carboxyl groups involved in 
the hydrogen bonding are approximately coplanar. 
The deviations of the carboxyl atoms from the least- 
squares plane through the guanidyl group (equation 
8) are 


O(1)—0-28, O(2)+0-06, and C(3)—0-05 A. 


Hydrogen bonding between arginine molecules also 
occurs parallel to the 6 axis in which O(2) from the 
carboxyl group is bonded to N(12’)H(14’) of the 
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yguanidy! group and N(9’)H(10') of the guanidyl group 
is bonded to N(5) of the amino group. The two 
hydrogen atoms, H(2) and H(3), in the amino group 
are not involved in hydrogen bond formation. This is 
a rather unusual occurrence. It is expected that 
molecules in a crystal usually arrange themselves to 
form a maximum number of hydrogen bonds (see, 
e.g., Fuller, 1959). 

The water molecules form infinite chains perpendic- 
ular to the arginine chains and parallel to the a axis 
by hydrogen bonding to each other. In addition, the 
water molecules form lateral hydrogen bonds to the 
arginine molecules. They are the donors of two 
hydrogen bonds and the acceptor for one. 

The hydrogen bond lengths and angles are listed in 
Table 5. The values are consistent with those found in 
many other crystals (Fuller, 1959). The four hydrogen 
bonds around W(1), the oxygen atom in H2Q(1), 
are arranged in a roughly tetrahedral configuration, 
whereas the three hydrogen bonds around W(2) are 
approximately in a plane 

The intermolecular approaches between 
non-bonded atoms are those between N(5) and 


N(12’) (3-42 A), and between N(5) and N(11”) (3-36 A), 


closest 
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Table 5. Hydrogen-bonds lengths and angles 
Bond 


N(9)H(10) 


Angle 


N(5’) 2-888 A C(8)N(9) 


C(10)N(9) 


N(5’) 112-7 
-> N(5’) 110-; 
C(10)N(12)--- O(2’) 135-6 
C(10)N(12)---O(1”’) 116-: 
C(1LO)N(11)--- O(2”) 121-: 


N(11)- Wi(1)- 
N(11)--- W(1)-- 
N(11)--- Wl): 
W(2)--- W(1)- 
W(2)--- Wl) 
W(2’)--- WL) 
O(1)--- W(2)--- 
O(1)--- W(2) 
Wil) W(2) 


N(12)H(14)--- O(2’) 2-897 
N(12)H(13) -- - O(1) 2-869 
N(LI)H(12) - - - O( 2”) 2-887 


W(1)H(15)---O(1’) 2-739 
N(LI)H(11)--- WL) 2-912 
W(2)H(18)--- W(1) 2-954 
W(1)H(16)--- W(2’) 2-709 
W(2)H(17)---O¢l) 2-890 


~W(2) 101-7 
O11’) 109-9 
W(2’) 138-0 
-OUl') 105-6 
> W(2’) 
-O(1') 
W(1) 
- Wil’) 
Wil’) 


where hydrogen bonds may have been expected to 
form, and 3-45 A between two water molecules not 
hydrogen-bonded to each other. 


We wish to acknowledge the invaluable assistance of 
Mr. Peter Gum, who performed all the high-speed 
machine calculations and programmed many of the 
routines for the IBM 7090 and IBM 7030 (STRETCH). 
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ABSTRACT Structures in four different crystal forms of 
{Leu'|zervamicin (zervamicin Z-L, Ac-Leu-Ile-Gin-Iva-Ile*- 
Thr-Aib-Leu-Aib-Hyp'®-Gin-Aib-Hyp-Aib-Pro'’-Phol, where 
Iva is isovaline, Aib is a-amino isobutyric acid, Hyp is 4-hy- 
droxyproline, and Phol is phenylalaninol), a membrane channel- 
forming polypeptide from Emericellopsis salmosynnemata, have 
been determined by x-ray diffraction. The helical structure is 
amphiphilic with all the polar moieties on the convex side of the 
bent helix. Helices are bent at Hyp’® from ~30° to ~45° in the 
different crystal forms. In all crystal forms, the peptide helices 
aggregate in a similar fashion to form water channels that are 
interrupted by hydrogen bonds between N*H(Gin'') and 
O*(Hyp") of adjacent helices. The Gin"' side chain is folded in 
an unusual fashion in order to close the channel. Space is 
available for an extended conformation for Gin'', in which case 
the channel would be open, suggesting a gating mechanism for 
cation transport. Structural details are presented for one crystal 
form derived from methanol/water solution: Cg<H49N ;g027 
10H,0, space group P2,, a = 23.068(6) A, b = 9.162(3) A, c = 
26.727(9) A, B = 108.69(2)° (standard deviation of last digit is 
given in parentheses); overall agreement factor R = 10.1% for 
5322 observed reflections {|F,| > 30(F)]; resolution, 0.93 A. 





Several acyclic, a-aminoisobutyric acid (Aib)-containing 
polypeptides of fungal origin (Fig. 1) form voltage-gated 
channels in phospholipid bilayer membranes (1-6). The best 
studied members of this class are alamethicins, 20-residue 
peptides for which a crystal structure has been reported at 
1.5-A resolution (3). Several models for alamethicin ion 
channels have been proposed (3-6), based on the helical 
conformation observed in crystals. The molecular packing in 
alamethicin crystals, however, did not directly provide a 
model for a membrane-phase channel aggregate (3). The 
zervamicins are a related class of peptides isolated from 
cultures of Emericellopsis salmosynnemata (7). These pep- 
tides are shorter (16 residues) than alamethicin and contain 
several polar residues (Thr, Gln, Hyp) distributed evenly 
throughout the sequence (8). A preliminary report suggests 
that zervamicins exhibit lower pore-forming ability and en- 
hanced antibacterial properties compared with alamethicin 
(8). More recently, the zervamicins have been shown to form 
voltage-dependent multilevel ion channels in bilayer mem- 
branes (ref. 9; M. S. P. Sansom, personal communication). 
Although the conductive activity is qualitatively similar to 
that of alamethicin, the kinetics of switching between states 
is appreciably faster for zervamicins. The major polypeptide 
components of zervamicin (IIA and IIB) contain a Trp 
residue at position 1. During HPLC fractionation, a minor 





The publication costs of this article were defrayed in part by page charge 
payment. This article must therefore be hereby marked *‘advertisement™ 
in accordance with 18 U.S.C. §1734 solely to indicate this fact. 


Alamethicin | (11) Ac-Aib-Pro-Aib-Ala-Aib” Ala(Aib)-Gin-Aib-Val 


Aib! °_Gly-Leu-Aib-Pro-Val! 5-Aib-Aib-Glu-Gin-Phol2° 


Antiamoebin I (II) Ac-Phe-Aib-Aib-Aib-D-Iva?-Gly-Leu-Aib-Aib-Hyp! 


Gin-D-Iva-Hyp(Pro)-Aib-Pro!5-phol 


Zervamicin IIA (IIB) Ac-Trp-le-Gin-Aib(Iva) lle>-Thr-Aib-Leu-Aib 


Hyp! °-Gin-Aib-Hyp-Aib-Pro!5-phol 


Zervamicin Z-L Ac-Leu-lle-Gin-Iva-lle>-Thr-Aib-Leu-Aib 


Hyp! 9-Gin-Aib-Hyp-Aib-Pro!5-Phol 


Fic. 1. Sequences of some Aib-containing antibiotics that trans- 
port ions across membranes. Phol, phenylalaninol; Iva, isovaline 
(a-ethylalanine); Hyp, 4-hydroxyproline. 


component containing Leu at position 1, [Leu']zervamicin 
(zervamicin Z-L), has been isolated. The component has 
been crystallized and the structure determined in four differ- 
ent crystal forms to resolutions as high as 0.93 A and R = 
10.1% (Table 1). In each of the four different crystal forms, 
space groups P2, and P2,2)2), the peptide molecules asso- 
ciate to form similar discontinuous water channels, which 
appear to give further insights into the nature of an ion 
channel, a possible gating mechanism, and the dynamics of a 
bent helix molecule. 


EXPERIMENTAL PROCEDURES 


{Leu']Zervamicin was isolated from a heterogeneous zer- 
vamicin mixture by HPLC ona reverse-phase Cg column, as 
described (10). Crystals were grown by slow evaporation of 
methanol/water solution or by vapor diffusion of water into 
ethylene glycol/ethanol solutions containing about 3 mg of 
peptide in 0.4 ml of solvent. X-ray diffraction data were 
measured with a four-circle automated diffractometer with 
CuK, radiation on crystals under the conditions listed in 
Table 1 for a best set of data. Lowering the temperature for 
some of the crystal forms resulted in shattering the crystal. 
Rapid scan speeds, 7-15°/min, were used in order to collect 
most of the data before significant deterioration of the crystal. 
The structures were solved by a vector-search procedure (11) 
contained in the SHELXx84 package of programs (MicroVAX 
version of the SHELXTL system of programs; Siemens Instru- 
ments, Madison, WI). The search model was based on the 
backbone atoms of the known structure of a 16-residue apolar 
analog of zervamicin (12) and was followed by a partial- 
structure procedure (13) by which the atoms in the side chains 
were located. The detailed structure only of crystal A will be 
reported here. Full-matrix anisotropic least-squares refine- 





Abbreviations: Aib, a-aminoisobutyric acid (a-methylalanine); Hyp, 
4-hydroxyproline; Iva, isovaline (a-ethylalanine); Phol, phenylalan- 
inol. 
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Table 1. Crystal parameters for Ac-Leu-lle-Gln-Iva-lle-Thr-Aib-Leu-Aib-Hyp-Gln-Aib-Hyp-Aib- 


Pro-Phol, CgsHy40N 1s022.xH2O 





Parameter Crystal A** 


Crystal B* 


Crystal C Crystal D 





Space group P2, 
Cell dimensions 
a,A 23.068(6) 
b,A 9.162(3) 
c,A 26.727(9) 
B, degrees 108.69(2) 
Growth solvent CH,OH/H,O 
No. of reflections 
(\F.| > 301 F)} $322 
Crystal condition —45°C, dry 
Volume per 
asymmetric unit, A* 2676 


21.857(4) 
9.381(3) 
26.744(6) 
105.22(2) 90 90 
CH;OH/H,0 


20°C, dry 


P2,2)2) P2)2;2) 
10.337(2) 
28.3897) 
39.864(11) 


10.160(3) 
28.252(9) 
40.338(9) 
HOCH,CH,OH t 


4239 2251 
20°C, in capillary 20°C, in capillary 


2924 2895 





*R = SIF - 


F.|\/=\F.| for all data with |F,| > 30(F) = 10.1%; resolution, 0.93 A. Estimated SDs for 


backbone: ~0.023 A (bonds); ~1.3° (angles). Estimated SDs for side chains: ~0.03 A (bonds); ~1.7° 


(angles) 


‘Crystals A and B were present in the same crystallizing vial. 
*Crystal D was obtained by soaking crystal C in aqueous KCI (no K* in crystal D); approximate 
coordinates for the peptide are known, but least-squares refinement was not completed. 


ment was performed on the C, N, and O atoms in the peptide, 
after which H atoms were placed in idealized positions with 
C-H distances of 0.96 A and allowed to ride with the C or N 
atoms to which each was bonded for the final cycles of 
refinement. The thermal factors for the H atoms were fixed 
at U;,, = 0.125. Water molecules for a total of 11 sites (one 
pair of sites is mutually exclusive) were located in a series of 
difference maps and added stepwise to the refinement. The 
phenyl group in the Phol residue was constrained to be a 
regular hexagon with C-C distances of 1.395 A, although the 
position and orientation of the hexagon and the isotropic 
thermal values for each atom in the hexagon were subjected 
to least-squares refinement. Finally, the amide moiety in the 
Gin’ side chain, particularly the N*(3) atom (in all the crystal 
forms) was particularly mobile; therefore the bond lengths 
and bond angles in the -C(O)NH)> moiety were constrained to 
reasonable values while the orientation of the amide group 
and the isotropic thermal factors associated with the C°(3), 
O*(3), and N*(3) atoms were allowed to vary in the least- 
squares process. With the above constraints, and not quite a 
full sphere of data to a resolution of 0.93 A due to the onset 
of crystal deterioration, the final R factor is 10.1% for 5322 
data observed with |F,| > 30(F). Bond lengths and bond 
angles (estimated SD, approximately 0.02 A for bonds and 
1.2° for angles) do not show significant or systematic differ- 
ences from expected values.* 


RESULTS 


The Peptide Helix. Both alamethicin and zervamicin have 
a high proportion of Aib residues, which are strong helix 
formers (14-16). Zervamicin has polar moieties in residues 3, 
6, 10, 11, 13, and 16. As expected for a helical structure with 
polar residues spaced at each third or fourth position (17, 18), 
all of these polar moieties are on one side of the zervamicin 
helix (Fig. 2), thus giving the helix an amphiphilic character. 
Although Gin'! would appear to break the rule, the torsions 
about the C*—C* and C°’—C® bonds in Gin" are —51° and 
-55°, respectively, which fold the polar -C(O)NH, group 
back from the nonpolar face to the polar face of the helix. The 
unusual conformation and the probable mobility of the Gin"! 





+Supplementary material consisting of atomic coordinates for C, N, 
and O atoms, bond lengths and angles, anisotropic thermal factors, 
and coordinates for hydrogen atoms has been deposited with the 
Cambridge Structural Data Base, University Chemical Laboratory, 
Lensfield Road, Cambridge CB2 IEW, U.K. Lists of observed and 
calculated structure factors are available from I.L.K. and J.L.F.-A. 


34 Naval Research Reviews 


side chain strongly suggest the location of the gating mech- 
anism during ion transport (see below). Additionally, the 
polar face is enhanced by the exposed backbone carbonyl 
oxygens in residues 6, 7, and 10, which do not participate in 
intrahelical hydrogen bonds (Table 2) due to the presence of 
Hyp”® and Hyp”. These carbonyls are available for hydrogen 
bonding to the neighboring water molecules or another 
peptide molecule or for making transient ligands with cations 
that are being transported. 

Despite the presence of two Hyp residues and one Pro 
residue, the backbone is entirely helical. Torsional angles are 
listed in Table 3. An a-helix is formed from the amino 
terminus through the ninth residue. A B-ribbon containing 
Hyp’, Hyp”, and Pro® is twisted into a helix that extends 
to the carboxyl terminus. The helix is bent at Hyp’®, giving 
the peptide an overall banana shape. It is interesting that the 


HYDROPHOBIC RESIDUES 


Fic. 2. Conformation of [Leu']zervamicin drawn with the coor- 
dinates obtained from crystal A. The C* atoms are labeled 1-16. The 
side chain in Gin" (filled circles) is not extended to the hydrophobic 
face but folded backward to the polar face of the helix. 
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Table 2. Hydrogen bonds 





Peptide helix 


N---Oor C=0:---N 
O---O,A i angle, degrees Donor Acceptor oO 
Head-to-tail N(1) os)? 2.93 . 122 W(1) o076)§ 
N(2) 0(16)* 2.96 : 140 Wil) 010) 
Intramolecular N(3) O*(3) 3.07 2.32 W(2) OG} 
5— 1 Ni4) O10) 85 : 167 Wwi(2) Wil) 
N(5) O(1) 15 158 Wi(3) 0*(13) 
N(6) O(2) 25 ] 146 WG) Ww(2) 
Intramolecular O”(6) O(2) 79 . W(4A) Wwi3) 
5—1 N(7) OG) 88 03 155 W(4B) WG) 
N(8) O(4) 3.36 131 N*(3)8 W(4B) 
N(9) O15) 2.98 2 165 Wi5) WI(4A) 
0°10)? (7) 2.75 : WS) W(6) 
N*(11)8 O16) 3.09 W(6) O14) 
N*(11)8 0°(10) 2.91 : Wwi6)' 0*(3) 
B-Ribbon N(11) O(8) 2.88 . O16) W(6) 
N(12) Or9) 3.15 ; wi?) Wi5) 
Intermolecular 0°(13) o*(11)* 2.67 8: W(8) O(11) 
B-Ribbon N(14) 011) 3.01 2. Ww(s) w9)" 
N(16) O13) 2.82 : w(9) O(12) 
wg) Wi8) 
Wwi10) O15) 


Water structure 





N---Oor 


Type Donor Acceptor 





Intermolecular 


BSRRRSSB 


= 
x 


E388 


wWwWwWNN WN NNN SH WS & WW iV N \ 


~~ 
Uw 





*The H atoms were placed in idealized positions with the N—-H distance equal to 0.96 A. 

*Symmetry equivalents 1 + x, y, 1 + z to coordinates listed in supplementary material (see footnote at end of Experimental Procedures). 
*Symmetry equivalents 2 — x, —% + y, 1 — z. 

§Symmetry equivalents 2 — x 

‘Symmetry equivalents x, 1 + 
Symmetry equivalents x, —1 


~wt+y.i4—z. 


conformation of the backbone of a completely apolar syn- 
thetic analog of [Trp']zervamicin (12) is almost identical to 
that of the backbone in the natural [Leu']zervamicin. The 


configuration of the Iva residue at position 4 is R, confirming 
an earlier stereochemical assignment (8). 

Channel Formation. The peptide molecules in all four 
polymorphs have very similar helical conformations, intra- 
molecular hydrogen bonds, head-to-tail hydrogen bonds, and 


interpeptide hydrogen bonds. The differences in conforma- 
tion are manifested mainly by the bend of the helix, 36° for A, 
30° for B, and 45° for C (structure D is not sufficiently well 
refined at this time for us to make a measurement). In all the 
crystals, the polar faces associate in a common antiparallel 
fashion. Fig. 3 illustrates this association with a superposition 
of the envelopes of the structures, showing the peptide 
molecules A, B, and C. In the middle, the molecules in each 


Table 3. Torsion angles* (degrees) 





Residue 1) Ww 


Leu! —51* —45 

lle? —63 —45 164, 

Gin? —69 —40 —67 

Iva‘* —57 —54 ] 57 

Ile* —74 —35 172, 

Thr® —§59 —44 —64, 

Aib’ —59 —37 

Leu® —9] —24 ] —74 —75, 

Aib? —55 —39 

Hyp” —63 —16 —24 6 —35 19% 

Gin"! —89 ~9 —5] ~54 —44, 

Aib”’ —§51 —37 —177 

Hyp» -77 ~6 —168 -21 34 —33 20** 
Aib'* —§3 —44 —178 

Pro’ —80 —13 —170 —10 24 —29 22** 
Phol!® —115 +724 ~69 -—75, 101 








The torsion angles for rotation about the bonds of the peptide backbone (@, w, and w) and about bonds of the amino acid 
side chains (y") are described in ref. 19. 
*Estimated SDs, ~1.5° 
*C’(O),N(1),C%(1),C'(1). 
$C(3),C%(3),C(3),0%(3), +15°; C9(3),C%3),C°(3),N"(3), —139°. 
§N(6),C%(6),C(6),07(6), —64°; N(6),C7(6),C*(6), C%6), 179°. 
§€9(10),N(10),C*(10),C9(10), 3°; C2(10),C4(10),C%(10),0°(10), —81°. 
C(11),C%(11),C%(11),0°(11), —44°; C9(11),C%(11),C%11), N*(11), 135°. 
**C9(13),N(13),C%(13),C9(13), 1°; C(13),C9(13),C%(13),0%(13), —83°. 
*¥C9(15),N(15),C%(15),CF(15), —8°. 
4N(16),C*(16),C(16)(CH2OH),O(16)(CH OH). 
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HYDROPHOBIC FACE 


Fic. 3. Schematic diagram showing a superposition of the enve- 
lopes of [Leu']zervamicin helices in crystal A (solid line), crystal B 
(dotted line), and crystal C (dashed line). The polar face of the helices 
associates in an antiparallel fashion and forms a discontinuous water 
channel that contains a different number of water molecules in each 
crystal form. The channel is closed in the middle by interpeptide 
hydrogen bonds involving backbone carbonyls O(6) and O(7) and 
side chain moieties O°(H yp”) and N*(Gin"'). The size of the mouth 
of the channel is determined by the movement of O(Phol’®) and 
N‘(Gin*) in adjacent peptide molecules 


crystal are held together with interpeptide hydrogen bonds 
O(Thr®) -- - N*(Gin''), N*(Gin'') --- O°(Hyp"®), and 
O*(Hyp"”) - - - O(Aib’). The banana-shaped molecules spread 
apart from each other at either end by differing amounts, 
related to the bend in each helix. The resulting cavities are 
filled with varying amounts of water depending upon the size 
of the cavity. For crystal A there are eight water sites (one 
disordered) at each end. These are shown in detail in Fig. 4, 
along with the hydrogen bonding between water molecules 
and between peptide and water molecules. 

The mouth of the water cavity is constricted or expanded 
depending upon the bend in the helix. The O(Phol"®) --- 
N‘(Gin*) separation is 3.7 A in crystal B, 6.8 A in crystal A, 
and >8 A in crystal C. Despite the flexibility of the peptide, 
atom O(Phol'®) is relatively rigid since it participates in head- 
to-tail hydrogen bonding. The -C(O)NH, end of Gin’, how- 
ever, is quite mobile and N‘(3) makes hydrogen bonds mainly 
to disordered or transient water sites. The carbonyl O(3) 
makes an intrahelical hydrogen bond to the backbone N(3) 
moiety of the same residue in crystal A but not in crystal B. 

Gating. A continuous water channel is not formed between 
the polar surfaces of adjacent peptide molecules. A view of 
the middle segments of adjacent molecules, where interpep- 
tide hydrogen bonding takes place and closes the channel, is 
shown in Fig. 5. The view is perpendicular to that in Figs. 3 
and 4 and is directed into the helices. Three molecules appear 
to define the closed channel. The side chain of Gin" is not 
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Fic. 4. Water molecules (@) in crystal A of [Leu']zervamicin. 
The two peptide molecules are related by a twofold screw axis. The 
atoms in the side chain of Gin! are marked with stripes. Hydrogen 
bonds are indicated by dashed lines. 


extended as Gin’ of alamethicin (3) but is folded backward to 
form a hydrogen bond with O*(Hyp’’). The Hyp residue is 
quite rigid. It is conceivable that under the application of a 
potential, the Gin"! side chain can open the channel by 


i eT ® 
oO O—O 4 > Y 
\. HyP10 “~~ # 
Ri a> % GLN11 
SB, 
r wi " 
x = es 


zy » 


y ‘ : 
X+-~ .e—__Oo7¢ 
»~< THRE 
x ies 
As 


ot 


Fic. 5. A possible gating mechanism by Gln! and Thr® in 
channels formed by [Leu']zervamicin. The view of the closed 
channel is perpendicular to the orientation in Figs. 3 and 4. Only the 
middle residues of three adjacent peptide molecules are shown. The 
water molecule closest to the closed portion of the channel is W(1) 
(checkered). The carbonyl oxygen in Aib’ and the hydroxyl in Hyp’® 
(striped) are reijatively rigid moieties and may define a boundary of 
a hypothetical open channel. 
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rotating into the adjacent space at the right (Fig. 5). If a 
column of water is pushed through the channel in a single file 
by a cation (20), it may also be possible that the side chain in 
Thr® can rotate about the C*-—C* bond to further enlarge the 
channel. As shown in Fig. 4, there are many carbonyl, 
hydroxyl, and amino moieties along the channel to provide 
ligands for a cation that is being transported. 

Nonpolar Face of Helix. Helical peptides forming a channel 
across a membrane would have their nonpolar side facing the 
membrane. The molecular recognition and aggregation of 
apolar helical peptides are not particularly selective. Such 
helical peptides have been found to crystallize with the helix 
axes either in an all-parallel motif or in an antiparallel motif. 
Even peptides with the same sequence have either all-parallel 
or antiparallel aggregation in different crystal forms (21, 22). 
Similarly in [Leu']zervamicin, in crystals A and B (space 
group P2,), the nonpolar face of the helix associates in an 
antiparallel fashion, while in crystals C and D (space group 
P2,2,2,), the nonpolar face associates in a parallel fashion. 
The inference can be drawn that the nonpolar face is also 
fairly insensitive to the helix axis direction with respect to the 
membrane. 


DISCUSSION 


[{Leu']zervamicin is a bent helical molecule, 29 A long, with 
all the polar side chains on the convex side and the nonpolar 
side chains on the concave side. Alamethicin (3), also a bent 
helix, has a much smaller degree of bending and has fewer 
polar moieties on the convex face (Fig. 1). The helix in both 
peptides is flexible as shown by different amounts of bending 
in the three conformers in the same crystal of alamethicin and 
the four conformers in four different polymorphs of [Leu']- 
zervamicin. Although complex irregular water channels were 
observed in the alamethicin crystals, the polar residues are 
not directed toward the center of the channel as might be 
expected for a membrane ion channel (3). In each of the four 
crystal forms of [Leu']zervamicin, the polar faces of adjacent 
peptide molecules form the boundary of a water channel. 

The molecular recognition of the nonpolar side of 
{Leu']zervamicin is not very specific. Either parallel or 
antiparallel association takes place in the different crystal 
forms. In contrast, the polar face is very specific in forming 
an antiparallel association with an adjacent molecule and a 
hydrogen bond between Hyp” and Gin", among others. The 
Hyp” - - - Gin’! hydrogen bond interrupts the water channel. 
The remarkable observation is that the Gin" side chain has 
to reverse its direction, in a swinging-arm fashion, to ap- 
proach the Hyp” of an adjacent, antiparallel peptide. It may 
be indicative of the gating mechanism in ion transport under 
the application of a potential. The size of the water channel 
near its interruption is restricted to one water molecule or one 
cation. The implications from the crystal structure analyses 
are that the water channel is not an open channel—i.e., that 
there is no passive diffusion of cations. This point has been 
emphasized further by soaking crystals in KCI solution for 
various times at various concentrations. Short of the crystal 
cracking and disintegrating, the soaking produced crystals 
with somewhat different cell parameters, but with no K* 
penetration. 

Time-resolved fluorescence anisotropy measurements on 
helical bilayer-spanning 21-residuc peptides have shown that 
the amplitude of fast fluctuations are largest at the helix ends 
and decrease toward the center of the helix (23). The four 
crystal structures show an expansion or constriction of the 
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mouth of the channel, 0(16) - - - Gin’, that varies from 3.7 A 
to >8.0 A. The cation may be trapped in the channel by 
constriction of the mouth (24). 

One of the interesting implications of the crystal structure 
of [Leu']zervamicin is the role played by the mobility of the 
long side chain of Gin residues: the demonstrated mobility of 
Gin’ in controlling the size of the mouth of the water channel, 
and the implied mobility of Gin"! (by its unusual conforma- 
tion) in probable participation in the gating mechanism. For 
proteins, the role of Gin residues has been characterized as 
relatively indifferent and has no extreme properties or violent 
preferences (25). Perhaps one of the distinctive functions of 
this residue is in the regulation of ion-transport peptides. 
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Looking Back 


What is Matter? Never Mind. 
What is Mind? No Matter 
Punch XXIX, 1855 


Dr. Isabella Karle, today, working on a crystal 
structure analyses project in her laboratory at 
NAL. Dr. Isabella Karle came to NRL in 1946 
from the University of Michigan with her hus- 
band Dr. Jerome Karle; where she was teach- 
ng freshman chemistry to engineering students 


Dr. Karle in 1964, the year she received the Navy 
Superior Civilian Service Award for outstand- 
ing achievements in molecular spectroscopy 


Naval Research Reviews 


Dr. Isabella Karle about 1948 at the electron diffraction 
apparatus. Soon after coming to NRL, Jerome Karle 
designed the large and complicated apparatus which 
Isabella assembled, aligned and fine-tuned and used 
for experimentation that determined three-dimensional 
structure and vibrational motion in gaseous molecules 








Dr. Karle aligning a crystal on an X-ray camera in the late 1960s 


In 1980, RADM Baciocco, Chief of Naval Research, pre- Dr. Fred Saalfeld, Director of the Office of Naval Re- 

sented Dr. Karle the Robert Dexter Conrad Award, the search, presenting the Secretary of the Navy Distin- 

highest honor for scientific achievement in the Navy guished Achievement in Science Award in 1987 to Dr 
Karle 
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Dr Karle (second from right) on the occasion of receiving the Federal Women’s Award in 1973 


Dr. Karle (right) when she received the degree of Doc- 


tor 
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f Humane Letters, Georgetown University, in 1984 


Naval Research Reviews 


Dr. Isabella and Dr. Jerome Karle have collaborated on 
many projects during their 47 years at NRL. It was Dr 
Isabella Karle who bridged the gap between the theory 
and experiment by translating mathematical results into 
practical procedures for determining the three dimen- 
sional structure of molecules 











The Bower Award 
of 
The Franklin Institute 


The Franklin Institute Of Philadelphia was founded in 1824 to promote 
scientific inquiry and to recognize scientific achievement. The Institute's Com- 


mittee on Science and the Arts was founded in 1834 and serves as the Institute's 
science selection committee. 

Through the careful work of this committee of scientists and engineers, the 
Franklin Institute was among the first to recognize the outstanding contribu- 
tions of such individuals as Marie Curie and Albert Einstein. The Bower Award 
and Prize in Science was made possible by a $7.5 million bequest to the Franklin 
Institute by Henry Bower (1896-1988), a Philadelphia chemical manufacturer. 
The award, which consists of a gold medal and $250,000, is presented annually 
to a distinguished scientist of the world for outstanding work in the life or 
physical sciences. There is a rigorous process of nomination and review, involving 
national and international boards of advisors. The selection mission for this 
world-wide award is to find a person who reflects the scientific and humanitarian 
genius of Benjamin Franklin. 

Aside from the Bower Award for Science, the Franklin Institute gives annu- 
ally a Bower Award for Business, and 16 different medal awards. 
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